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[. INTRODUCTION 


A. FLOW SEPARATION 

Real fluids flowing over bluff bodies separate under the influence of 
adverse pressure gradient as momentum is consumed by wall shear. Separa- 
tion is best understood in terms of its consequences. In fact, it was The 
keen observations of separation in diffusers That led Prandtl to introduce 
the concept of transition or boundary layer theory. Even though the 
boundary layer theory has revolutionized fluid dynamics, the phenomenon 
which gave impetus to its discovery still remains an enigma. 

Through the years various characteristics of separated flow about 
bluff bodies in general and about circular cylinders in particular have 
been partially understood through observations, measurements, and analyses. 

The separation point may be mobile as on a circular cylinder or fixed 
along a line as on a wedge or blunt-based body. When a separation point 
is mobile, the pressure decreases up to the point of maximum velocity 
(slightly upstream of separation) and then ¡increases up to the point of 
separation. When the separation point is fixed, the pressure continues 
To decrease up to the point of separation. Consequently, there are signi- 
ficant differences in the separation of flow about a circular cylinder 
and a sharp-edged body. In either case the velocity distribution, rate 
of vorticity flux, and the amplitude and frequency of oscillation of the 
Separation points are not independent of the wake in the downstream side 
of the body. It is this interaction that makes the calculation of sepa- 
rated flows exceedingly complex. In other words, one needs the evolution 


of the boundary layer on the forebody (upstream of separation) of the 








cylinder in order to determine the instantaneous position of the separa- 
tion points. But, the characteristics of the boundary layer are, in turn, 
strongly affected by the fluid motion in the wake. Evidently, one needs 
some kind of an interactive mechanism whereby the effect of the wake on 
the boundary layer, and hence separation, can be given proper considera- 
tion. In mathematical terms one can state that the unseparated potential 
flow, or the d'Alembert flow, does not even constitute a first approxima- 
tion to the separated real flow. Partly the realization of OS raci, 
partly the difficulties associated with the integration of the Navier- 
Stokes equations across a singularity (separation point), and partly lack 
of sufficient computational capability have prevented the prediction of 
the kinematic and dynamic characteristics of separated flows about bluff 
bodies. In the following, the basic characteristics of flow about circu- 
lar oylinders are described with special emphasis on those aspects which 
prepare the reader for the understanding of the numerical mode! developed 
during the course of this work and the assessment of the similarities and 


dissimilarities between the predictions and observations and measurements. 


B. VORTEX SHEDDING 

The separation of flow gives rise to alternate vortex shedding. This 
Phenomenon is Intrinsic to the flow itself and has been known at least 
Since the times of Leonardo de Vinci. In 1878, almost exactly 100 years 
ago, Strouhal discovered that there is a relationship between the frequency 
Of vortex shedding, the velocity of flow, and the diameter of the cylinder. 
This relationship, which is now known as the Strouhal Number, St, is 
given by St = fD/U. Strouhal's experiments gave an average value of St = 
0.185. During the course of numerous experiments since then, it has been 


experimentally demonstrated that: (i) the Strouhal number is dependent on 





the Reynolds number and, as first suggested by Lord Rayleigh, must de 
written as f = St(Re)U/D; (ii) the vortex shedding frequency is only an 
average quantity for a given flow and that there are always secondary and 
tertiary frequencies present; and finally, (iii) there are regions of flow 
such as the region of drag crisis where it is quite difficult to identify 
a particular frequency. It is evident from these facts that in some re- 
gions of flow one must speak of the spectral content of the wake fluctua- 
tions in assessing the value of the Strouhal number and the forces exerted 
on the cylinder. The foregoing discussion of the Strouhal number raises 
two questions: firstly, is it possible to devise a universal Strouha! num- 
ber which would be applicable to all two and three-dimensional bluff 
bodies, and secondly, is it possible to theoretically predict the Strouha! 
number tnrough, for example, the perturbation of the d'Alembert flow. An 
attempt along these lines has been made by Sacksteder LI] who found St = 
0.2028 for an infinitely large Reynolds number. The reasoning behind 
This prediction and its validity will not be discussed here further. 

The former question of the possibility of devising a universal Strouhal 
number nas been of interest to various researchers. Table |, shown below, 


lists those definitions together with their originators. 


TASCE || 


Body Type 



























FAI /UK two-dimensional] bodies Roshko [2] 
EAU two-dimensional! bodies Bearman L3] 
SS axisymmetric bodies Calvert [4] 


f 2 4 U two-dimensional bodies Simmons [5] 





In Table |, h' represents the wake width as determined from Roshko's 


notched- hodograph theory [2]; &, lateral spacing of the vortices as 





determined by Kronauner's minimum drag criteria as reported by Bearman [3]; 
2 , the wake width ata distance of the wake establishment region; and K = 


vI-C where C 


pb is the base pressure coefficient. All these definitions 
require either the measurement or the solution of one or more characteris- 
tics of the wake. Thus, their capability to predict a vortex shedding fre- 
quency is quite limited. 

It is advantageous to note at an early stage those factors which in- 
fluence the Strouhal number. Firstly, small amplitude transverse oscilla- 
Tion of cylinders or cables distinctly regularize the vortex shedding and 
hence the Strouhal number. For example, the experiments conducted at The 
National Physical Laboratory [6] have shown that the Strouhal number even 
in the drag crisis region remains perfectly constant and equal to about 
0.21. Secondly, the vortices do not peel off from a cylinder as perfectly 
two-dimensional line vortices. As it was first shown by Humphreys [7], 
There is a cell pattern separated by a correlation length. The longer the 
cylinder the larger the number of cells. The consequences of the forma- 
tion of such cells are that there are phase shifts in pressure along the 
cylinder and hence a reduction in the integrated transverse force. In 
other words, the transverse force calculated through the integration of 
The pressure distribution at a given section is not identical to the aver- 
age transverse force determined experimentally. All numerical analyses 
which model the complex three-dimensional flow as a two-dimensional flow 
must necessarily predict a transverse force or lift coefficient somewhat 
larger than those found by measurements. Other consequences of the cel! 
formation concern the strength of the vortices. The vortices are not 
necessarily oriented in a plane normal to the flow and the measurements 


(flow visualization or electro-optical techniques) yield only that 





component of vorticity which is normal to The plane of flow. Additional- 
ly, experimental results necessarily depend on the length-to-diameter 
ratio, end conditions, intensity and scale of turbulence in The ambient 
flow, the rigidity of the cylinder mounting, and of course on Reynolds 
number and the surface conditions of the cylinder. Most recent experi- 
ments (Shih and Hove [8]) have shown that roughness affects the base 
Pressure in the same manner as the small transverse oscillations (beyond 
a threshold limit of about 0.10 diameters). Uniform roughness yields a 


perfectly uniform base pressure and removes the cell pattern. 


ERESI STANCE IN STEADY FLOWS 
|. In-Line Forces 
The understanding of the fluid forces exerted on bluff bodies 

comes primarily from measurements. Of all bluff bodies, the circular cy- 
linder has been subjected to the largest number of investigations. Those 
experiments which constitute milestones in the determination of the magni- 
tude and characteristics of resistance are the works by Wieselsberger [9], 
Fage and Johansen LIO], and Rhosko LII]. A plot of the drag coefficient, 


D 
Reynolds number defined by Re = UD/v. It is a well known fact that for 


C~ defined by Co = 2F/oLDU“, is shown in Fig. | as a function of the 


Reynolds numbers larger than about 107 and smaller than about 2x10?, the 
drag coefficient for a circular cylinder remains essentially constant at 

FE Sluse of about 1.2. For Re > 2x 107, the shear layers emanating from 

the separation points become turbulent. Since a turbulent flow can ex- 
change larger momentum between the inner and outer flow (i.e., between the 
boundary layer and external flow), it can sustain a larger adverse pressure 
gradient and even reattach to the cylinder. The reattachment of the tur- 


bulent shear layers gives rise to the formation of so called separation 
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bubbles and to the displacement of the attached flow further downstream 
from the first separation point. The reattached flow cannot continue To 
remain attached indefinitely and reseparates at an angle of about |30 de- 
grees from the front stagnation point. Thus, the transition of the flow 

in the shear layers from a laminar to a turbulent state results in a con- 
Siderably smaller wake and hence in a considerably reduced drag coefficient. 
This particular flow regime is commonly called the drag crisis region or 
the critical region. As the Reynolds number increases, parts of the bound- 
ary layer adjacent to the first separation point become turbulent and the 
second separation point begins to move gradually upstream. This region is 
called the transcritical region. At higher Reynolds numbers, the entire 
boundary layer becomes turbulent and the drag coefficient reaches a nearly 
constant value of about 0.6. There is considerable scatter in the data 

in the supercritical region because of the difficulties associated with 
flows at such high speeds. 

The present knowledge of the drag coefficient for a smooth cylin- 
der is limited to Reynolds numbers smaller than about 10’. There is great 
practical interest in the determination of force coefficients, i.e., drag 
coefficient, lift coefficient, and the Strouhal number for projects such 
aS ocean thermal energy conversion (OTEC), large cooling towers, high rise 
Structures, etc. A careful consideration of the various possible methods 
To obtain data at Reynolds numbers in the order of 108 or ja shows that 
There are practically insurmountable difficulties. The existing wind and 
water tunnels are not capable of yielding such high Reynolds numbers with- 
out the merce effects of blockage and compressibility. The pop-up and 
drop-down experiments with large cylinders are constrained by cavitation, 


weight, and stability conditions. It might be possible in the future to 


ZI 








close the high Reynolds number gap by building large cryogenic wind tun- 
nels operating with freon or liquid helium. 

The drag coefficient depicted in Fig. | may be influenced by a 
number of factors such as roughness and turbulence. Numerous investiga- 
tors, notably Fage and Warsap [12], Achenbach [13], and Guven [14], have 
clearly demonstrated that the larger the surface roughness, the sooner is 
The occurrence of the critical region and the higher the dip in the drag 
coefficient (see Fig. 2). Turbulence plays a similar role in the evolu- 
tion of the drag coefficients. I+ should be noted from Fig. 2 that the 
drag coefficient for a rough cylinder at supercritical Reynolds numbers is 
about twice that for a smooth cylinder. Furthermore, the experiments by 
Güven [14] and Szechenyi [15] have shown that at sufficiently high 
Reynolds numbers there is a Reynolds number independent region. This ex- 
perimental fact has been termed the independence principle and used (Shih 
and Hove [8]) to extrapolate the drag coefficient data to Reynolds numbers 
in the order of 10 If the validity of the independence principle can 
be demonstrated unambiguously, it may obviate the need for difficult ex- 
periments in cryogenic tunnels. 

2. Transverse Forces 

Lift or transverse~force coefficient data are far less numerous and 
exhibit a great deal of scatter. Figure 3, which represents most of the 
existing lift coefficient data, shows that the lift coefficient, C, defined 


L 
by C = (amplitude of lift SS O ÉS may vary from 0.2 to 1.0 in the 
subcritical Reynolds number range. In the critical range it is impossible 
to define a lift coefficient without resorting to its root-mean-square (rms) 


values. It is for this reason that some of the data appearing in Fig. 3 


eee in terms of the rms values of the lift coefficients. For a more 
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detailed discussion of the spectral content of lift in the critical region 
Ee Bublitz (17). 

The reasons associated with the large variations of the lift co- 
efficient are directly related to the mobility of the separation points, 
formation of vortex cells, lack of spanwise correlation, end effects, tur- 
bulence present in the ambient flow, etc. Since each experiment is carried 
out with a different length-to-diameter ratio, blockage ratio, and end 
conditions, the number of vortex cells over the cylinder and the effect of 
the flow between the end of the cylinder and the Tunnel wall give rise to 
widely different lift coefficients. 

The effect of roughness on the lift coefficient is to reduce its 
scatter and to increase its magnitude. As noted earlier, roughness dis- 
tinctly regularizes the vortex formation and the base pressure thereby 
eliminating the effect of phase shifts which cause the said scatter in 
lift. Data supporting the above experimental facts wil! not be reported 
here for sake of brevity. For additional details the reader is referred 


to the work by Szechenyi LI5]. 


D. RESISTANCE IN TIME-DEPENDENT FLOWS 

There are an infinite number of time dependent flows and it is impos- 
sible to discuss their characteristics in any unifying manner. One might 
Obtain a clearer insight into such flows by considering only those cases 
which are relatively more manageable for numerical or laboratory experiments. 
Among such flows which have been subjected to extensive numerical analysis 
and experiments are the impulsively started flows, uniformly accelerating 
Flows, and harmonically oscillating flows. Obviously, no distinction is 
made here between the said fluid flows about a body at rest and a body 


undergoing similar time dependent motions in a fluid otherwise at rest. 
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|. Impulsively Started Flows 


The impulsively started flow has been by far the one most investi- 
gated numerically through the use of the Navier-Stokes equations. These 
investigations (Collins and Dennis [18], and Yang [19]) have been confined 
to the very early stages of the motion during which the separation points 
move from the rear stagnation point to about 109 degrees. An exception 
to the foregoing is the work reported by Thoman and Szewczyk [20]. They 
have solved the Navier-Stokes equations through the use of a suitable 
finite difference technique for various Reynolds numbers and found that 
the drag coefficient overshoots for a relative flow displacement (normalized 
time) of Ut/c = 5 and then gradually reaches steady state values. 

Only two separate investigations of the impulsively started flow 
have been reported. The first is that carried out by Schwabe [21] upon 
Prandtl's suggestion and the other by Sarpkaya £22, 23]. Schwabe's experi- 
ments were conducted with a 9 cm diameter cylinder at a Reynolds number 
Of about 600. His drag coefficient, calculated indirectly through the use 
of motion pictures, rapidly reaches a value of about 2.0 and remains con- 
stant thereafter. Schwabe's findings are not in conformity with either 
the numerical calculation or with the experimental results obtained by 
Sarpkaya. Figure 4 shows three sets of representative data at Reynolds 
numbers from 16000 to 35000 for the drag coefficient. Clearly, the very 
initial instants of the motion cannot be analyzed because of the problems 
associated with the generation of an impulsively started flow. At Ut/c 
Slightly greater than zero, the force acting on the cylinder is almost 
entirely inertial. As the separation sets in at the rear stagnation point 
and rapidly moves upstream, the drag also increases rapidly and Co reaches 
a value of about 1.5 at Ut/c between 4.0 and 5.0. Then Cy drops rapidly 


To the steady state value for smooth cylinders in steady uniform flows at 
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Ut/c larger than about 7.0. Detailed flow visualization experiments by 
Sarpkaya [22] have shown that the reason for the overshoot in Ch is the 
Symmetric development of two vortices in the cylinder wake and the rapid 
accumulation of vorticity in each vortex. If must be noted that for such 
a symmetric vortex configuration, the only means by which The vorticity 
may be cancelled are tne exchange of vorticity between the two vortices 
and the cancellation of vorticity with oppositely-signed vorticity, so 
called counter-vorticity, generated in the rear boundary layers. In other 
words, the rate of accumulation of vorticity in the vortices is consider- 
ably larger than that at any other time. Other vorticity cancellation 
mechanisms which come into play during the separation of the vortex from 
its feeding sheet, as proposed by Gerrard [24], do not play any role dur- 
ing the growth of the symmetric vortices. Thus, the rapid accumulation 

of vorticity gives rise to the drag overshoot. As the asymmetry sets in, 
vortices shed alternatingly and the flow gradually becomes asymptotically 
Steady. The possible mechanisms which cause the asymmetry in an otherwise 
symmetrical flow about a symmetrical body are discussed in greater detail 
later. Suffice it to note that one of the fundamental facts of fluid 
Tlows is that symmetric causes do not give rise to symmetric effects, as 
most eloquently philosophized by Birkhoff [25]. 

Following the onset of asymmetry, the cylinder experiences a 
Transverse force, as shown in Fig. 5, in terms of the lift coefficient and 
The normalized time Ut/c. The lift force fluctuates with the Strouhal fre- 
quency. The shedding of each vortex gives rise to a smal! fluctuation in 
Che The frequency of this fluctuation is twice the Strouhal frequency 
Since two vortices (one from top and from bottom) are shed during a given 


Eycle of lift oscillations. 
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2. Other Time-Dependent Flows 


As noted earlier, the uniformly accelerating flows and harmonical- 
ly oscillating flows have been subjected to numerous analytical and experi- 
mental investigations. Among those notable are the investigations carried 
out by Iverson and Balent [26], Keim L27], Hamilton [28], Odar and Hami | ton 
[29], Sarpkaya and Garrison. [30], and Sarpkaya [31]. The basic question 
has been and continues to be the formulation of a resistance equation which 
takes into consideration the effects of acceleration, velocity, and the 
history of the motion. All E have shown that resistance in 
time-dependent flows in general and in harmonically oscillating flows in 
particular cannot be considered as a juxtaposition of resistance in ET 
flows at corresponding instantaneous Reynolds numbers. Stokes, in his 
celebrated paper on the sinusoidal oscillation of a sphere in a fluid 
otherwise at rest, has shown that the fluid resistance is comprised of a 
linear velocity-dependent force and an acceleration-dependent inertial 
force. The coefficients for both forces depend on the Reynolds number and 
The frequency parameter D7/VT where T is The period of oscillation. Ob- 
viously, Stokes! analysis is valid for only very small, as yet unspecified, 
Reynolds numbers. Inspired partly by Stokes! contribution and partly by 
Practical considerations, Morison [32] proposed a resistance equation which 


may be written as 


a dU 
F = 5 0A, Cp Ulul + o¥ C, == (1) 


in which U represents the instantaneous velocity; Are the projected area; 
¥, the displaced volume; and Cp and Cu’ the drag and inertia coefficients 
respectively. For harmonically oscillating flow about a stationary cylin- 


der where the fluid velocity is U = "Uy cos 8 (here 8 = 2nt/T) Eq. (1) 
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reduces To 


A 
| zè 


TT ; 
0,170 Cm sin 8 CZ) 


= - Cp [cos 8|cos 8 + 
The essence of this equation is that the resistance in harmonic flow is 
assumed to be equal to a linear combination of a velocity-squared dependent 
drag force and an acceleration-dependent inertia force. In the midst of 
insufficiently clear insight into the kinematics of the complex flow, if 
has not yet been possible fo demonstrate The reasons for the differences 
between the measured force and that calculated from Morison's equation for 
Uy T/D values in the vicinity of 19. Additional experiments and numerical 
analysis, possibly through the use of the discrete vortex model, might shed 
further light on the subject and ¡into the overall understanding of resis- 
tance in Time-dependent flows. It is evident from the discussion of both 
the steady and unsteady flows that analytical or numerical methods are 
urgently needed with which relatively Inexpensive investigations can be 
conducted to clarify the complex wake-boundary-layer interaction mechanism. 
Of all the possible flows which might be considered as a candidate for such 
an analysis, the impulsively started flow is certainly the one most manage- 


able and most fundamental to the understanding of ail other flows. 


See METHODS OF NUMERICAL ANALYSIS 

With the advent of high speed computers, numerical methods began to 
Play a greater role in fluid dynamics. Currently four major methods of 
analysis are in use: finite difference, finite element, marker-cell tech- 
nique, and the discrete vortex method. 

Finite difference methods have been used for the solution of flow 


about a circular cylinder for relatively small Reynolds numbers (Re = 100) 
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either through the use of the Navier-Stokes equations as they are commonly 
written or through the more frequently used vorticity equations. Ingham 
[33] calculated the steady flow past a circular cylinder at Re = 100. 
Thoman and Szewczyk [20], as noted earlier, carried out a similar numeri- 
cal analysis of flow over stationary and rotating cylinders at Reynolds 
numbers as large as oe. There are a number of such examples in the 
literature but They are not cited here because of their lack of relevance 
to the subject under consideration. 

The finite element method, which came into prominence first in the 
area of structural! analysis and Then permeated through heat transfer and 
fluid mechanics, has been primarily used for unseparated flows about 
bluff bodies. As far as fluid mechanics is concerned, the outstanding ex- 
amples of such applications may be found in two recent volumes of work 
dealing with finite elements in fluid mechanics [34]. - Bretanow and Ecer 
[35] made an attempt to calculate the separated flow about a circular 
cylinder at smal! Reynolds numbers but the results have been inconclusive. 
Applications of the finite element method to flows about foils where the 
flow remains attached have yielded results in conformity with experiments 
(Bretanow and Ecer [35]). 

The marker-cel | technique developed by the researchers at the Los 
Alamos Scientific Laboratory during the past 20 years uses the finite 
difference form of the Navier-Stokes equations and then adjusts the pres- 
sure in each cell to insure that the equation of continuity is satisfied. 
Numerous ingenious techniques have been developed by Fromm and Harlow 
[36] and Nichols and Hirt [37] for the efficient adjustment of pressure 
in each cell. Such techniques will not be dealt with here. The researchers 
Of Los Alamos Scientific Laboratory have applied the marker-cell technique 


To the flow past cylinders, spheres, rectangular blocks and to various other 
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free-surface flows. These applications dealt with very low Reynolds num- 
bers (in the order of 100). Even though the technique has considerably 
increased the understanding of the kinematics of the wake, the precise 
nature of the time-dependent force acting on any one of the bodies cited 
above remains obscure. 

The discrete vortex model, which is included as one of the four major 
models, is the subject of extensive discussion in the next chapter. There- 
fore, the remainder of this section is confined to a discussion of the suc- 
cesses and shortcomings of the first three methods. 

All numerical models are by their very nature approximate. They satis- 
fy certain critical conditions only at discrete time steps. Secondly, they 
require mesh or element sizes which vary from one region of flow to another 
to insure that the characteristic gradients of the flow (such as velocity 
or pressure) do not significantly or abruptly change over any one element. 
For example, the finite difference method when applied to the viscous flow 
Over a circular cylinder requires that the size of the elements on and near 
The cylinder be considerably smaller than the boundary layer thickness. 
OTherwise, the evolution of the boundary layer, the positions of the sepa- 
ration points, er evolution of the wake cannot be accurately calcu- 
lated. Since the boundary layer thickness for laminar flow is proportional 
to 2 one realizes that only for relatively small Reynolds numbers can 
a sufficiently large boundary layer thickness and a correspondingly large 
mesh size be obtained. This fundamental difficulty has prevented the 
finite difference, finite element, and the marker-cel! technique from making 
Significant progress into the calculation of steady or time dependent flows 
about bluff bodies at Reynolds numbers of practical interest. None of the 
Three techniques cited above has been applied to the solution of flow prob- 


lems associated with oscillating bodies. The reasons for this are rather 
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self-evident since the solution requires not only a moving grid but also 
a grid whose mesh ie must be varied with time to minimize the effect of 
sharp velocity and pressure gradients. Finally, it should be noted that 
all three numerical techniques present "stiffer" solutions and introduce 
an artificial numerical viscosity whose influence cannot be readily 
assessed. Not to be entirely critical of the said methods, for there are 
corresponding shortcomings of the discrete vortex model, one must note 
That in the wide spectrum of Reyno!ds numbers there will always be some 
interest in low Reynolds number flows. Thus, such numerical techniques 
will continue to play some role in separatéd flows where the discrete 
vortex mode! is not advantageously suited. It therefore appears that no 
one model wil! stand out as superior and universally applicable and that 


all methods will continue to complement each other. 
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Il. FUNDAMENTALS OF THE DISCRETE VORTEX MODEL 


A. INTRODUCTION 

The discrete vortex model (DVM) is a potential flow representation of 
the shear layers emanating from the separation points. Any separated vis- 
cous flow is completely specified if the spatial and temporal distribution 
of vorticity in the entire flow field is known. For reasons described in 
the previous chapter, it has not yet been possible to achieve this objec- 
tive through any exact or approximate methods. IT has been observed and 
experimentally verified by Fage and Johansen LIO] that the shear layers 
at sufficiently high Reynolds numbers are quite thin and the vorticity is 
confined into regions of tightly spiralled vortex sheets. This fact alone 
invites one to divide the shear layer into a number of small segments and 
concentrate the vorticity of each segment into a line vortex situated at 
the geometric center of the segment. This idealization permits one to 
use the powerful complex function theory to determine the kinematics and 
dynamics of the flow. Evidently, the validity of the results can only 
be judged in light of experimental facts. 

In the following, only the gross features of the method are described 
and the details are presented following an historical résumé. 

At a given time the separated shear layers are represented by vortices 
situated at the center of each segment. The potential flow theory permits 
one to calculate the velocity of each vortex induced at its center by all 
Other vortices and the ambient flow. Then the vortices are convected for 
a short time interval A+ through the use of a suitable convection scheme 


TO new positions. It is then necessary to introduce new vorticity into 
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the flow at or near the separation points to account for the vorticity 
generated in the boundary layer. The introduction of new vortices (nascent 
yortices) represents an addition to each set of vortex arrays. Then the 
process is repeated. 

In spite of its apparent simplicity, there are a number of problems 
which must be resolved in carrying out the said calculations: (i) the vor- 
ticity introduced at each time step must be related to the characteristics 
of the boundary layers upstream of the separation points; (ii) one must 
insure that the point vortices continue to represent the vorticity at each 
segment in accordance with the fundamentals of hydrodynamics; (iii) allow- 
ance must be made for the interaction of oppositely-signed vorticity and 
the decay of vortices; (iv) methods must be devised to limit the number of 
the point vortices to minimize the computation time without affecting the 
Kinematics and dynamics of the flow; and (v) procedures must be established 
for the flow to retain its inherent flexibility without introducing 
artificial stiffening effects (e.g., to allow for the mobility of stagna- 
tion and separation points). 

The foregoing conveys the idea that the discrete vortex model is and 
will continue to rely on a strong interaction between potential and vis- 
cous flow theories. Neither is capable of providing a solution alone but 
Together they can make inroads towards describing the basic features of 
the separated flow provided that sufficient ingenuity is incorporated into 
the said interaction. Lastly, it must be noted that al! the numerical 
methods described above depend on the availability of a high-speed digital 


computer. 


B. HISTORICAL DEVELOPMENT 
The DVM has been in existence for the past 46 years since its incep- 
Tion by Rosenhead [38] in 1932. Rosenhead investigated the so-called 
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"Helmholtz instability" of a surface of discontinuity, separating two 
fluids of equal speed and density but moving in opposite directions. His 
primary objective was to give The shear layer between the two fluids an 
initial disturbance by situating the shear layer along the curve y = 0.| 
sin (21x/A) and then observe the evolution of the sheet with time. To ac- 
complish this, Rosenhead represented the shear layer by point vortices 
placed along the sheet at intervals such that their projections along the 
x-axis were equi-distant. Each vortex was assigned a strength of UA/6n 
where n is meer of vortices per wavelength,A, and U the speed of 

the two fluids. Keeping the vortex strengths constant, Rosenhead followed 


the deformation of the sheet by convecting the vortices with a simple 


ep een ed 
t+At T 


Bi the j-th vortex; q, the computed velocity vector; and At, the time 


Eulerian integration of z quat where pa is the position vector 
interval. Using n = 2, 4, 8, and I2, Rosenhead observed that in each case 
the sheet rolled up smoothly into concentrated clusters of vortices posi- 
tioned at intervals of A along the x-axis. Rosenhead's results, which were 
computed using a desktop calculator, are shown in Fig. 6. The tendency of 
a distribution of point vortices to roll up in a manner consistent with 
observations is clearly evident. 

Rosenhead's results remained undisputed until 1959 when the computa- 
tional power of the digital computer was applied to the analysis. Birkhoff 
and Fisher [39] repeated Rosenhead's analysis using a Runge-Kutta integra- 
tion scheme, 22 vortices per wavelength, and smaller time intervals. In- 
stead of improved results they observed that the vortices took irregular 
Paths and no longer rolled up smoothly. Birkhoff and Fisher concluded that, 
given sufficient time, an originally organized array of equi-strength point 
vortices would eventually become random and that the validity of represent- 


ing a shear layer by point vortices is questionable. 
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Shortly thereafter, Hama and Burke [40] undertook a similar computer 
study and verified the observations of Birkhoff and Fisher. However, Hama 
and Burke noted a conceptual error in Rosenhead's original method of plac- 
ing the discrete vortices along the shear layer. Specifically, if the 
initial vortex sheet was fo have a constant vorticity per unit distance 
along its length and be represented by equi-strength vortices, then the 
vortices should be placed at equal distances measured along the sheet vice 
at equal x increments. Using this arrangement Hama and Burke repeated 
Their calculations and found that the onset of irregular vortex motion was 
Significantly delayed. Eventually, however, a random pattern did result. 
Thus, at This point it was hypothesized that the onset of instability 
could be delayed although the questions of why it occurred and how it 
could be controlled were still a mystery. 

Without fully realizing the consequences of sheet instability (the 
result of which is discussed later), the DVM continued in its evolution 
relying on the assumption that meaningful results could be obtained up to 
the onset of irregular vortex motions. In 1962 Abernathy and Kronauer [41] 
applied the method to investigate the motion of two shear layers of oppo- 
sitely-signed vorticity separated by a distance ny and each given an 
initial perturbation as in Rosenhead's work. They observed that the vor- 
tices accumulated into clusters resembling the pattern of a vortex street. 
An example of their results is shown in Fig. 7. 

Encouraged by the demonstrated ability of the DVM to represent a vor- 
tex street, the next logical extension of the method was to analyze the 
flow past a bluff body. However, the inclusion of a solid body in the 
flow field necessitated new considerations. Specifically, the previous 
analyses had begun with a discretized sheet having a preassigned initial 


position and vorticity distribution along its length (usually taken as 


De 








Fig. 7. Example of Abernathy 
and Kronauer's calculations 
(AA ORSA nA]. 
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constant). The application of the DVM to flow past a bluff body required 
that the discretized shear layers be generated continuously by adding 
additional point vortices into the flow at regular time intervals and at 
locations representative of where the shear layers actually emanate from 
the body. The requirement for some method to account for vorticity genera- 
tion and introduction in an otherwise potential flow model posed a signi- 
ficant problem. Compounding the difficulty was the time dependent nature 
of both the generation rate and its point of introduction. 

The first DVM analysis to include a bluff body was performed by 
Gerrard [42] for a circular cylinder. Although a circular cylinder was of 
great practical interest it was also the most difficult one to analyze 
because of the unsteady vorticity generation and mobile separation points. 
Gerrard was forced to adopt a rather arbitrary method of introducing new 
vortices which severely influenced the results of his analysis, although 
it still represented a noteworthy application of the DVM. Noting the dif- 
ficulties associated with the mobile separation points and the a 
introduction, Sarpkaya [43] introduced a mobile separation region repre- 
sented by an infinite number of vortex sheets of vanishingly small vorticity 
connected to the nascent vortex. This analysis was restricted to ¡mpul- 
Sively started symmetric flow about a circular cylinder. In the ensuing 
years many investigators bypassed the difficulty introduced by the mobility 
of the separation points and dealt with bodies with sharp edges (e.g., 
inclined plates, blunt-based bodies, wedges, etc.) where the separation 
Points are fixed. Despite the growing number of applications of the DVM, 
a number of disturbing difficulties remained unresolved and a full inter- 
action between the boundary layer, separation points, and the wake had 
not been established. These difficulties required the use of a number of 


nOndisposable parameters. Consequently, the proper balance between 
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heuristic reasoning and computational formulation had to evolve through 
the recognition of the reasons underlying the basic difficulties. 

Briefly, the most important of these difficulties stemmed from the 
proximity of the Wen eas 46 each other and to the solid boundary, and 
from the necessity of representing the continuous process of vorticity 
generation at the singular separation STE by a discrete process. The 
occasional proximity of two point vortices resulted in large mutually- 
induced velocities and hence the shooting of the two vortices rapidly 
away from the flow field. This was particularly true when the vortices 
ceased to be located at the geometric center of the sheet segment which 
they were supposed to represent. Furthermore, the spiralling of the sheets 
brought vortices on two adjacent spirals close to each other and gave rise 
to a traveling and growing instability. Thus, it was not possible to 
correctly represent the core of Tightly spiralled sheets. The proximity 
of The vortices to the boundary and hence to their images inside the 
body resulted not only in similar instabilities but also in unrealistic 
local pressures. This, in turn, resulted in unnatural fluctuations in 
aag and [ift forces. 

Faced with the problems just cited and the desire to enhance the 
applicability of the method, various investigators devised several numeri- 
cal cures. Chorin and Bernard [44] introduced a cut-off distance from 
The center of each vortex beyond which the vortices were to behave like 
potential vortices. For radial distances smaller than the cut-off distance 
the vortices were assumed to exhibit a solid body rotation. Their calcu- 
lations did not, however, result in an improvement of the analysis of 
separated flow about a circular cylinder. Moore [45], in dealing with 
Two vortex sheets emanating from a plane lifting surface, introduced the 


method of amalgamation. When a vortex sheet was spiralled tightly enough, 
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the last vortex on The sheet was connected to a core vortex. At each Time 
step the last vortex was amalgamated into the core vortex increasing its 
strength. This procedure partly avoided the difficulties stemming from 
the proximity problems of the adjacent coils of the spiralling sheet. 

Both investigations dealt with the devising of methods to overcome the 
instability problems rather than with the understanding of the real cause 
Ehe difficulty. 

The basic question of why an initially equidistant set of vortices 
should after several convections come in close proximity of each other was 
resolved by Fink and Soh L46] who demonstrated that the primary source of 
the difficulty was in failing to rediscretize the sheet at each time in- 
terval of the calculation. In other words, it is necessary to insure that 
the point vortices are always located at the geometric center of the seg- 
ment of the sheet which they are supposed to represent. The complex con- 
jugate velocity alz.) ata point a due to a single vortex of strength T 


K 


located at a iS given by. 


T 
qlz,) kan 
er 


(5) 


For a series of N discrete vortices on a vortex sheet, the same velocity 


at a point Z becomes [46] 
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KE ematical details of the relevant potential flow theory will be 
Presented in Chp. IIl. Here no attempt is made to introduce the precise 


‘definitions of the variables appearing. 


43 





Equation (4) resembles Eq. (3) for the induced velocity field due to point 
vortices except for the logarithmic term. The point a lies within the 
segment SS 541 /2? SAA without necessarily bisecting it at all times, and 
S measures the distance along the sheet. The consequences of Eq. (4) are 
that: (i) if the equivalent point vortex is not placed at the midpoint of 
its segment through rediscretization of the sheet at each time interval, 
then the logarithmic term does not vanish and the computational error in- 
creases depending on the number of vortices, the time interval used, and 
the total time of computation; (ii) the vortices, which initially bisect 
the segment which they are to represent, do not continue to do so at the 
Succeeding time intervals; (iii) The use of finite vortex cores, amalgama- 
Tion of vortices at the center of the spiral, or other techniques only 
delay or minimize the accumulation of the errors resulting from the loga- 
rithmic term in an amount related to the distance between El and the center 
of the segment; and that (iv) the growth of the computational error may be 
significantly reduced by placing each discrete vortex at the midpoint of 
its segment, i.e., by placing the vortex at q = ea at each 
Time interval. Only through such a procedure can one make the logarithmic 
term vanish. 

The calculations are then carried out at each time step by represent- 
ing the vorticity density by an entirely new set of equidistant vortices 
whose strengths are adjusted to give a good representation of that density. 

The foregoing procedure does not resolve all of the computational 
errors particularly in regions where the radius of curvature of the sheet 
is small. Fink and Soh applied this technique to the early stages of 
relatively simple flows with fixed separation points. Their examples did 


not deal with the alternate shedding of vortices. 
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Even though the method of rediscretization essentially resolved the 
problems stemming from the proximity effects, two fundamental links between 
the potential flow model and the behavior of real fluids remained unre- 
solved. The first of these is the relationship between the vorticity to 
be introduced at each time step and the boundary layer which generates it. 
The second is the decay of vortices by various mechanisms. The first 
question was initially considered by Deffenbaugh and Marshall L47] who in- 
vestigated the impulsively started flow about a circular cylinder (without 
recourse to rediscretization). They used Pohlhausen's boundary-layer ap- 
proximation method to determine the position of the separation point and 
The vorticity to be introduced at each time interval. Their analysis pre- 
dicted a separation angle of SE = 64 degrees, which is considerably smaller 
Than the commonly accepted experimental value of a. = 8l degrees. - 

The question concerning the decay of vortices has remained practically 
obscure. Although several researchers employed methods to cancel opposite- 
ly-signed vortices whenever they came within a prespecified distance of 
each other, no systematic investigation of the effect of decay of vorticity 
was conducted. Often recourse was made to Prandt!i's suggestion that indi- 
vidual vortices ina vortex street retain only about 60 percent of the 
vorticity generated in the boundary layer during a period required to shed 
a single vortex. Some investigators, in efforts to obtain results in 
partial conformity with Prandtl's suggestion, automatically applied a 40 
percent reduction fo the strength of each nascent vortex. Seldom was men- 
Tion made of the experimentally observed fact that the rate of decay of 
vortex strengths varies with distance behind the cylinder. 

The foregoing represents the state-of-the-art prior to the beginning 
of the present investigation. A relatively more detailed account of the 
historical evolution of the DVM is found in Clements and Maul] | [48] and 


Fink and Soh [46]. 
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Ill. THEORETICAL AND PHENOMONOLOGICAL FOUNDATIONS 
ORRE DISCRETE VORTEX MODEL 

A. INTRODUCTION 

This chapter deals with the characteristics of inviscid and viscous 
flows about a circular cylinder from the forward stagnation point to the 
vortices in the far wake. Mathematical principles, boundary layer con- 
cepts, and relevant hypotheses are introduced in preparation for the dis- 
cussion of the details of the discrete vortex model developed in the 


present study. 


B. FLOW KINEMATICS 


The complex potential for a two-dimensional inviscid flow is given by 
w(z) = @ + i 


where 6 and Y represent, respectively, the velocity potential and the 
stream function. The velocity components which satisfy the Cauchy-Riemman 


conditions are given by 


EBEN er a ea 
u a Sy and v Sy ae (9) 


In terms of the complex potential, one has 


dw(z) E 


FE - u + iv (6) 


The flow about a circular cylinder is generated by combining a uniform 
flow with a doublet at the center of the cylinder (see e.g., Mi Ine-Thomson 
[49]). This yields 


w(z) =w= =- Ulz + T (7) 
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in which -Uz represents the uniform flow and the remainder the doublet. 
The introduction of a vortex into the flow field exterior to The cylinder 
requires, by virtue of the circle theorem, the introduction of an image 
vortex interior to the circle. Thus, the flow about a circular cylinder 


with N number of real vortices is given by 


N 


2 
Sf un. o Y ae 
w= -Ulz + =) + DE r (&2n(z a A = )) (8) 


u n 
n= | 


It should be noted that a real vortex located at z = Z has an image at 
z = A and at z = 0. The image at the center of the circle, which 
gives rise to terms like P2n(z), must be excluded due to the condition 
that the vortices shed from the cylinder leave a circulation opposite to 
their own on the cylinder. Clearly, the absence or presence of an image 
vortex at the center of the circle does not alter the boundary condition 
(dQ/an| = 0) since the streamlines of a vortex are circles. No re- 
z|=c 
striction has been imposed regarding the time dependence of the velocity 
U other than that it be uniform. 
Substituting Eq. (8) info Eq. (6), one obtains the velocity af any 


point z as 


N 
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Equation (9) may also be used to evaluate the convective velocity of a vor- 
tex at z = Zo provided that the singular term T /(z-z,) is omitted by virtue 


of the condition that a vortex cannot impart a velocity on itself. 
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It is advantageous to work with dimensionless parameters by introduc- 


ing the following change of variables: 


IE a iv! = (uFfiv)/U, TT" = PUC (10) 


w!(z) = w(z)/Uc nie 


Substituting into Eqs. (8) and (9), letting U = 1.0, c = 1.0, and dropping 


the primes for sake of simplicity, one has 
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where u, v, T and z are now the normalized variables and w the normalized 
complex potential function. Hereafter, al! measures of distance, time, 
velocity and circulation are given in terms of the above normal ized 
variables. Any exceptions are clearly specified. 

Equation (12) is used extensively in convecting the vortices and in 
Calculating the lift and drag forces. In fact, about 90 percent of the 
computation time of the present model is consumed in the evaluation of 
complex velocities. This points out the necessity of keeping the number 
of discrete vortices as small as possible while maintaining a sufficiently 
small time interval to insure that the characteristics of flow are faith- 
fully reproduced. The selection of the time interval, like the selection 


Of the mesh size, is one of the subjective items of the method and 
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requires a careful parametric analysis. This is addressed to in a later 
section. 

In general, the calculation of velocities and the convection of vor- 
tices are straightforward and noniterative. Under certain circumstances 
vortex-vortex proximity between adjacent spirals and the vortex-boundary 
proximity may require special attention. The resolution of such secondary 
issues is deferred to another section since they are not of any major con- 


sequence either theoretically or practically. 


See FORMULATION OF THE RESISTANCE EQUATIONS 

The forces acting on a stationary circular cylinder of unit length by 
a time-dependent flow containing N vortices of strengths I and locations 
z can be calculated from the normalized form of the generalized Blasius 


theorem (Milne-Thomson [49]) to give 


Cyt IC, => INTA A AE 


where Co and C, represent the drag and lift coefficients respectively. 
Inserting Eq. (II) into Eq. (13) and evaluating the residues, Sarpkaya 
[50] has shown that 
i i 
+ — — - 
a > P tiv -v ) = Cu vl (14) 
n= | 

in which u and ve represent the velocity components of the n-th real vor- 
Tex and u and Vo, the velocity components of the image of the n-th vortex. 


Equation (14) may be written in a more compact and computational ly 


more suitable form as 


SR 
CSE DE > A (15) 
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where the dot denotes derivative with respect to time such that 
z =q =u + iv (16a) 
and 


nee y (1 6b) 





: q 
AE 
an = 2 (17a) 
Z 
n 
or by 
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The normalized force given by Eqs. (14) or (15) can also be evaluated in- 
directly through the use of surface pressures and the unsteady form of 
Bernoulli's equation. In either case the results are identical and the 
use of one method or the other is simply one of personal choice and pro- 


gramming efficiency. 


D. GENERATION AND INTRODUCTION OF VORTICITY 
The vorticity per unit length of a Thin shear layer separating two 


streams with velocities U, and U, is equal to the difference of the velo- 


cities U =U. This vorticity is transported with the average velocity of 
The two streams, i.e., CU, + U,)/2. This is a continuous process and the 


amount of vorticity transported during a time At becomes 


AT = (0-0, (0, +U l 


Ä ) At/2 = CU, 


2 
2 =U, Paak (18) 


50 





The discrete vortex model represents the continuous generation and flux 

Mi vorticity into the outer flow by introducing point vortices at discrete 
time intervals. The strength and position of introduction of each vortex 
must depend on the interaction of the unsteady vorticity distribution in 
the wake and on the appropriate boundary-layer considerations necessary to 
determine where the separation occurs. Those applications of the DVM 

which consider the separation points fixed Ignore not only the said inter- 
action but also artificially constrain the evolution of the flow. 

In the present analysis the instantaneous position of the separation 
points is determined through the use of an approximate boundary layer 
analysis. At a given instant the velocity distribution on the forebody 
of the cylinder is calculated and the corresponding pressure distribution 
is assumed To be:impressed on the boundary layer. This fundamental feature 
Of the boundary layer theory together with the integral momentum equation 
of von Karman permits the use of a number of separation prediction methods. 
Among the most notable of such methods are those given by Pohlhausen [51], 
RR ford [52], Thwaites [53], Curle and Skan [54], Timman [55], and 
Loitsianski [56]. All of these methods exhibit varying degrees of sensi- 
tivity to the external velocity distribution and may not al! yield identical 
or uniformly consistent results. 

It is important to note that the combination of potential and viscous 
flow methods cannot establish a one-to-one relationship between the Reynolds 
number and the flow characteristics. It is only the kind of separation 
criteria used (laminar or turbulent) that determines whether the flow is 
Subcritical, critical, or turbulent. The fact that the use of the DWM be- 
comes increasingly more accurate with decreasing shear layer thickness 
makes one realize that the results obtained with a laminar separation 


Criteria correspond to moderate to high subcritical Reynolds numbers. 
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The integral-momentum methods such as that of Pohlhausen require that 
the velocity outside the boundary layer be given and that the separation 
points be either fixed or fluctuate with small amplitudes and velocities 
about a mean position. In other words, the changes in the outer flow and 
the separation points should be such that the flow may be treated as 
quasi-steady. 

The outer velocity distribution may be calculated either numerical ly 
Through the use of an interactive scheme between the wake and the boundary 
layer or through the use of experimentally determined pressure distribution. 
Hiemenz [57] used the measured pressure distribution about a circular 
cylinder and approximated the velocity distribution at the edge of the 


boundary layer by 
> 3 2 
MU SEE A E O AS (19) 


where € represents the arc length measured along the cylinder from the 
Stagnation point. Hiemenz's solution indicates that separation occurs at 
an angle a. = 82 degrees, measured from the stagnation point. This cele- 
brated result has often been taken as proof that in subcritical flows 
Separation occurs at SØ = 82 degrees. The fact has not been realized that 
The entire velocity distribution and in particular that portion between 
The maximum velocity and the separation point (i.e., from 9 = 70 to 9 = 

80 degrees) dramatically affect the position of the separation point. In 
fact, Meksyn [58] has shown that the higher order derivatives of Hiemenz's 
velocity profile do not satisfy certain fundamental criteria regarding 

The separation point. At present there is no completely satisfactory velo- 
city distribution for the forebody of a circular cylinder either in sub- 
critical or in the supercritical regimes. Additional discussion of this 


Subject is found in Meksyn Gea) 
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Numerous experiments have yielded laminar separation points ranging 
from 8. = 70 degrees to 8. = 90 degrees. The commonly accepted value is 
about 80 degrees. Fage and Johansen's LIO] experiments, and others since 
then, have shown that the separation point is in reality a separation zone 
which extends over a finite region and that spanwise and chordwise fluctua- 
tion of This region do not permit a precise determination of the mean 
position, amplitude, and period of oscillation of the separation point. 

The purpose of the foregoing brief discussion is fo point out both the 
theoretical and experimental difficulties encountered in the determination 
of the separation point and the inadequacies of the appropriate laminar 
boundary layer methods. The prediction of the separation of turbulent 
boundary layers is even more uncertain (Stratford [52], Takada [59]) and 
are not discussed herein. {nstead, attention is focused on Ponlhausen's 
method for time-dependent boundary layers. The former is more appropriate 
to the flow conditions where the motion of the separation point and the 
velocity at the edge of the boundary layer may be treated as quasi-steady. 
Schun's [60] method is more appropriate for the initial instants of motion 
where the separation points move rapidly from a = |80 degrees to about 
9 = IIO degrees. The motion of the separation point commences after a 


5 


relative fluid displacement of (Schlichting L61_]) 


a | : 
Be ee a (20) 


The reason for this is that in impulsively started flows about rounded 
bluff bodies the separation does not begin until the fluid moves some dis- 
Tance past the body. The only exception to that is the flow about bodies 
with sharp corners in which case the separation begins at t = 0.0. 

For the case of impulsively started flow about a cylinder, the time 


in the boundary layer calculations is measured from the instant of the 
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inception of separation. Thus, a relative time of t = 0.35| must be added 
to that used in the boundary layer calculations in order to match the 
starting times. 

Pohlhausen's method is well described in the literature (Schlichting 
[61]) and only the highlights of the method are presented herein. Assuming 
a fourth order polynomial for the velocity distribution within the boundary 
layer, Pohlhausen has shown that von Karman's integra! momentum equation 


may be reduced To 
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Other variables appearing in Eqs. (21), (22) and (25) are 


UCE): the normalized velocity on the cylinder 


&: the arc length measured along the cylinder from the 
stagnation point 


ô: the boundary layer thickness normalized by the cylinder 
radius 


V: viscosity normalized by the ambient velocity and the 
cylinder radius (sometimes referred to as the unit 
Reynolds number). 
The values of K and F(K) are tabulated in terms of A (Schlichting [6l]). 
Similar tables are provided for other boundary layer approximations cited 
Previously. 


Integration of the coupled Eqs. (21) begins at the front stagnation 


Point and proceeds along both sides of the cylinder independently until 


54 





the separation criteria (A = -12) is satisfied for each layer. This re- 
sults in two distinct separation points which in general are not at 
identical angular distances from the mobile stagnation point. 

As noted earlier, the foregoing method is applicable only to the rela- 
tively established flow where the velocity at the edge of the boundary 
layer and the separation point do not significantly vary with time. Ouring 
the early stages of motion, however, the time dependence of the said para- 
meters is very strong and must be treated through the use of the time 
dependent Integral boundary-layer equations. Schuh [60] has extended 
Pohlhausen's method for the general calculation of the unsteady boundary 
layer. The details of his method are not presented here. The results 
obtained by Schuh for an impulsively started flow about a circular cylin- 
der are shown in Fig. 8. These results, based on the classical potential 
flow velocity distribution about a circle, show that the separation point 
rapidly approaches a value of about a. = |09 degrees. Schuh's method is 
not capable of incorporating the growth of the symmetric wake and the 
subsequent motion of the separation points and thus is not valid at large 
times. Thoman and Szewczyk [20] employed the finite difference form of 
The full Navier-Stokes equations and obtained results identical to those 
Shown in Fig. 8 for the initial instants of flow. They have further shown 
That the separation point once having reached an angular position of about 
| IO degrees moves gradually to the final value of the separation angle 
(about 80 degrees). Partly to confirm this conclusion and partly to assess 
the significance of time dependence on the instantaneous position of the 
Separation points, two sets of calculations are performed. Both use the 
time-dependent velocity distribution available through the model and both 
begin with an initial separation point of a. = 109 degrees. The first cal- 


culation uses Pohthausen's steady-state separation criteria with the in- 
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Fig. 8. Separation angle versus time [60]. 
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stantaneous velocity distribution. The second calculation uses the unsteady 
integral-momentum equations given by Schuh [60], including those terms 
associated with the unsteady velocity distribution. The results obtained 
from both sets of calculations are shown in Fig. 9. Evidently the two cal- 
culations yield almost identical results, pointing out the fact that the 
separation point moves rather slowly and that the contribution of the time 
dependent terms in the unsteady form of The equations are negligible. It 
is because of this reason that all calculations are based on Pohlhausen's 
separation criteria assuming the flow to be quasi-steady. This assumption 
improves as one approaches steady state. Figure 9 indicates that the 
angular speed of the separation point is about |0 degrees per second. The 
Same speed about the mean position of the separation point is about 2 de- 
grees per second. (As is shown later, the separation point oscillates 
about a mean position of a. = 7/7 degrees with an amplitude of 3 degrees 
and a period of about 10.) Based on the foregoing, an amplitude of oscil- 
lation as much as 15 degrees can be performed without the inclusion of 
Time-dependent terms. 

In the present model the calculations begin at the time of t = |.0 
when the separation point is at 8. = |09 degrees. It should be noted 
That a straightforward application of Pohlhausen's method with a potential 
flow velocity distribution (U = 2 sin 8) yields a separation angle of 
Eo = 109.5 degrees. 

Once the position of the separation points is determined, the rate of 


generation of vorticity is easily determined by writing 


6 
<= Sy 


37 





110 6; 


N sa = Une ready CompUTar on 


100 ——— Steady computation 





Fig. 9. Comparison of steady and unsteady 
senatet calculations. 
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where 6 is the boundary layer thickness; u, the velocity within the bound- 
ary Ber; and y, the distance along the outward normal. The contribution 
of dv/dx is excluded from the vorticity expression in conformity with the 
approximations of the boundary layer theory. 

The integration of Eq. (24) to obtain the circulation flux at a separa- 
tion point yields 
YO” 
man 


(25) 


aja 
15 


where Ue is the tangential velocity at the predicted separation point. 


Ordinarily the rate of change of circulation should be written as 


die 
ee OS 625) 


where u and U, represent respectively the velocities at the outer and 
inner edges of the separating shear layer. Extensive experiments with 


various bluff bodies by Fage and Johansen LIO] have shown that U, is about 
3 percent of u and that Sr is correctly represented by y A2. This is a 


universally accepted experimental and theoretical fact. 


In numerical calculations ar is replaced by AT/At and the strength of 


The nascent vortex is calculated from 


U Ent 
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The placement of nascent vortices into the flow to adequately represent 
The separating shear layers has been a source of difficulty. This is pri- 
marily because the discretization of a continuous process near a singular 
region (separation point) is an exceedingly complex problem. In the present 
analysis, no new procedures are needed for the determination of the strength 
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and position of the nascent vortex because the separation point is obtained 
from the boundary layer analysis. By virtue of its definition, the 
separation point must be the point at which both the velocity and The 

shear stress on the cylinder are q Thus, the nascent vortex of known 
strength must be placed on a radial line passing through the separation 
point at a distance m, i.e., at z , = (Ite) exp Li(m-6.)J, to render the 
tangential velocity on the cylinder at the separation point equal to zero. 
This yields a unique value for € as 


rl 
E A | (28) 


a zi 


The use of the boundary layer theory and the separation condition uniquely 
specifies the strength and position of the nascent vortices and allows the 
separation points to move along The cylinder as required by the interaction 
between the wake and the boundary layer. IT is of importance To note That 
the no-slip condition is satisfied only at discrete Time steps. Efforts 
to satisfy the said condition in a more continuous manner are counterproduc- 
tive because the decrease in At calls for unacceptably large computation 
Times. One of the essential features of the model is the bracketing of a 
Suitable interval for At, which provides an adequate scale of discretiza- 
tion without excessive computation time. 

The method described above properly accounts for the vorticity generated 
On the forebody of the cylinder. However, the flow in the wake of the 


Cylinder develops a boundary layer on the afterbody which generates 


En 


GÅ lengthy discussion may be made regarding the difference between the 
definitions of steady and time dependent separation criteria. However, a 
numerical analysis such as the DVM does not warrant the consideration of 
The distinction between the results stemming from the use of different 
separation criteria. 
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countervorticity. Procedures by which this relatively small amount of 
vorticity generation may be accounted for are far from simple. The velo- 
city distribution on the afterbody has strong spatial and temporal gra- 
dients due to the immediate proximity of the fluctuating wake. These 
Gradients render useless the separation prediction methods described 
ear | ier. 

Davis [62] attempted to account for countervorticity effects by using 
four separation points (two for the forebody and two for the afterbody 
boundary layers) but his method required the assumption of fixed separation 
points. Deffenbaugh and Marshal! [47] approximated the rear shear layer 
effect by assuming that the boundary layers remain attached over the same 
fraction of a region of adverse pressure gradient as that of the forward 
boundary layers. However, in view of the complex velocity distribution 
On the afterbody, it is incorrect to apply the same boundary-layer con- 
Siderations as were used on the forebody. It is safe to say that the 
numerical perturbations introduced into the model! in attempting to account 
for the secondary effects of countervorticity may cause more error than 
would the neglect of the effect altogether. This consideration plus the 
introduction of the circulation reduction mechanism (yet to be described) 
obviate the need for the separate accounting of the effect of the rear 


Shear layers. 


E. CONVECTION AND REDISCRETI ZATION MECHANISM 
The vortices emanating from a separation point may be convected in a 
number of ways. Three most promising schemes are those given by the 


following expressions: 


z(t+At) = z2(t-At) + 2q(t) At (29a) 
2(teAt) = z(t) + of3q(t)-q(t-At)] At (29b) 
z(t+At) = z(t) + q(t) At (290) 
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where q(t) is the complex velocity of the vortex evaluated at the positions 
and times indicated. These schemes yield more or less the same accuracy 
and the use of one in lieu of the others is primarily governed by the sim- 
plicity of the method and computation time considerations. In the present 
analysis the simple Eulerian scheme given by Eq. (29c) is used because of 
its simplicity and because the use of relatively smal! time intervals (AT = 
0.125) make the method sufficiently accurate. Clearly, the results obtained 
with the three schemes cited above, or with others, cannot be compared 
without taking into consideration the effect of the time intervals used 

and the discretization process (if one is employed). All such convection 
Schemes must assure, within the overall limitation of the model, that the 
vortex paths in regions of very small radius of curvature do not signifi- 
cantly deviate from those calculated with much smaller time Intervals. 

It has been explained in connection with the discussion of the histori- 
cal development of the DVM that all investigators, with the exception of 
Fink and Soh [46], let the vortices move in the form of vortex clouds and 
Took special measures to eliminate the difficulties arising from vortex- 
vortex and vortex-cylinder proximity effects. Not fo be entirely critical, 
The use of vortex clouds can yield force-transfer coefficients nearly 
identical to those obtained with the use of rediscretized sheets. However, 
One drawback of the vortex clouds is that it does not allow the formation 
Of spiralling vortex sheets as often observed experimentally. There are 
Other features of the rediscretization which will become clear later. 

The vortex clouds stem partly from the accumulation of the logarithmic 

errors and do not In any way represent a Turbulent vortex. In the early 
Stages of the present investigation the vortex clouds were used to obtain 
some idea about the gross features of the asymmetric wake. Figure 10 is 


a sample computer plot of vortex clusters. It is not possible to identify 
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Fig. 10. Flow configuration without rediscretization. 


A 


Pigs vies, flew configuration with rediscréetization. 
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a spiralling sheet in Fig. 10, but the shedding of the vortices and the 
thinning of the connecting vortex sheets between two adjacent vortex 
clusters are clearly observable. 

The fact that the vortices must represent each segment of the shear 
layer as accurately as possible and thus avoid the buildup of logarithmic 
errors in their convection gave rise to the method of rediscretization or 
redistribution. A sample flow pattern obtained with rediscretization is 
shown in Fig. II where the point vortices are connected by line segments 
to denote the position of the sheet. The mathematica! foundations of the 
concept have already been introduced. In the following, the details of 
The mechanism where the rediscretization is affected are described. 

Consider a segment of the vortex sheet at an instant t and assume that 
the vortices are connected by straight line segments (see Fig. 12). The 


distance between any two vortices is given by 





Fig. 12. Vortex sheet before rediscretization. 


o JZ, - NN (30) 


The length of the sheet from the position Z, (core of the spiral) to the 


n-th vortex is given by 
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n 


n 
s=) sò ET (31) 
j=2 


j=2 


The total length of the sheet for N number of vortices becomes 
Sy = > = (32) 


The vorticity per unit length of the sheet at each 2, is given by 


T 
_ n 
Veen E 
nti n-| 
Ne FÒN (33) 
Sa l',/òs, 


The circulation per unit length, Y,» may be tabulated as a function of SI 
through the use of Eqs. (31) and (33). This tabulation is then used as a 
piecewise continuous function denoted by y = y(s). The function y(s) uses 


a linear interpolation of the tabulated values, i.e. 


YOY 
ee so «s«s (34) 
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y(s) represents circulation per unit length and is an invariant at a given 
Time of the rediscretization process. 

The next step in the analysis is the replacement of the existing un- 
evenly spaced N vortices with N evenly spaced new vortices. In genera | N 
may be different from N. Replacement of N vortices is accomplished by 


starting at the core of the spiral, Zis and placing vortices at intervals 


Of ds measured along the straightline segments connecting the original 
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vortex positions (see Fig. 15). It should be noted that when The method is 
applied to the example depicted in Fig. 13, Zs and 2. are positioned such 
that ES - Z| = ås and ES > lè: la - Z,| = 68. In other words, dis- 


tances are measured along the straightline segments of the existing sheet. 





Fig. 13. Vortex sheet after rediscretization. 


The strengths of the new vortices are calculated through the use of the 
distribution function y(s) such that the circulation per unit length at a 
new position, a is The same as for the original sheet at the said posi- 


tion, l.e., 


I = y(o) ôs, 
A S 1705, _,| A 
Bo - SAA ER [MSN es) 
A = Y (Sy? ds 
where 
s* = (n=l) és 


The final step in the basic procedure is to insure that there has been no 


Net gain or loss of circulation during the rediscretization process. The 
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circulation deficit is computed as 


N N 
AT = >> r- yo i (36) 
n=] = | Ä 


n 


and an equal portion of the difference is added algebraically to each vor- 


mex, i.e., 


a 


r =f peL Erp 
n n a 
N 
Two parameters yet to be specified are ós and N. The present model 
does not change the number of vortices during rediscretization and thus 


N =N. 8S and Ñ are related by 


A 


Ss = S/N zn (38) 


It is important to note that N is the total number of vortices in a detached 
Sheet (The definition of a detached sheet is ¡introduced later.) and The 
rediscretization process does not change the position of either the first 


Or The N-Th vortex. For the vortex sheets attached fo a nascent vortex, 


\ 
b 


however, N represents the number of vortices excluding the nascent vortex. 
The reason for this is that the conditions (strength and location of the 
nascent vortex) prevailing at or near the separation point should not be 
altered by any process once the nascent vortex has been placed in the flow 
in accordance with the no-slip condition. Thus, for an attached sheet the 
rediscretization process stops at the vortex immediately adjacent to the 
Beiscent vortex. 
It is appropriate af this time to give some ¡dea about the actual 


Magnitude of some of the parameters mentioned above: 
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(1) the strength of a nascent vortex = 0.05 to 0.15 
(ii) N = 3 to 80 depending on the length of the sheet 
ei) és = 0. 
The above numbers are representative of those encountered for a stationary 


cylinder with a time step, AT = 0.125. 


F. VORTEX SHEDDING AND CIRCULATION REDUCTION 
I. Vortex Shedding 

One of the least understood of the wake formation processes is vor- 
tex shedding where the sheet connecting the spiralling vortex to the sepa- 
ration point is cut. That portion of the sheet which remains attached to 
the cylinder begins to roll up into a new vortex spiral while the detached 
portion moves downstream to become part of the vortex street. The shedding 
process is accompanied by a number of interactive phenomena such as the 
motion of the shed vortices, adjustment of the longitudinal and lateral 
Spacing of the vortex cores, and the reduction of circulation through the 
exchange of ly signed vorticity. These phenomena are so inter- 
locked that attempts to explain the reasons for the occurrence of one be- 
come in turn questions to be explained. Evidently one possible way to 
resolve this dilemma is to have an exact solution of the problem. Numeri- 
cal models such as the finite difference technique, having failed to 
elucidate the interactive mechanisms, left researchers with the need to 
develop techniques which could be used in more phenomonological models 
such as the DVM. It is absolutely essential that one begin the analysis 
with a model and then examine its validity through its predictions. 

Inspired by flow visualization and wake measurements, Gerrard [63] 
Proposed That the growing vortex spiral draws the opposite shear layer 


across the wake when it becomes sufficiently strong. The approach of 
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oppositely signed vorticity cuts off further supply of circulation to the 
vortex. The approximate time of the cutting of the sheet ¡s then considered 
as the beginning of the shedding process. This mechanism is sketched in 
Fig. 14. 

The wake comprised of such shed vortices continuously adjusts it- 
self in both lateral and transverse spacing. In general it is agreed that 
the wake may be divided into three parts. The first is the so called for- 
mation region Re = 5); the second, the stable region where the vortices 
exhibit the characteristics of a fairly uniform vortex street (5 Soke ES 
12); and the third, the unstable region extending beyond x larger than 
about 12 [64]. 

The lateral spacing of the vortex cores gradually ¡increases within 
the stable region and, like the longitudinal spacing, depends on the 
Reynolds number. 

Even Though considerable experiments have been conducted (particu- 
larly at low Reynolds numbers) to delineate the basic features of the wake, 
very few attempts have been made to measure directly the strength of the 
vortices as a function of x at a given Reynolds number. Admittedly, this 
is an extremely difficult undertaking. Several investigators (Schaefer 
and Eskinazi L64], Griffin [65], and Davies [66]) measured the longitudinal 
and transverse spacing of the vortices in the stable region and calculated 
indirectly the strength of the vortices (assumed to be identical) through 
The use of von Karman's vortex street model. Consequently, the vortex 
Strengths are approximate and highly dependent on the initial inputs. The 
most important conclusion resulting from such investigations is that the 
strength of the ee: decreases with increasing distance downstream. 

The question concerning how much of the circulation generated remains in 


a vortex has been the subject of discussion since the time of Prandtl. As 
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noted earlier Prandtl suggested that the net circulation remaining in a 
vortex is about 50 percent of That generated. Discrete vortex models used 
by Clements [67] and Sarpkaya L68] give a reduction of about 15 percent. 
This value is Significantly ler than the estimates of 40 percent re- 
duction (Fage and Johansen [10]), 57 percent (Roshko LI! ]) and 70 percent 
(Bloor and Gerrard [69]). I+ should not be inferred from the foregoing 
that the A area hans acquired a net circulation of about 40 per- 
cent continue to retain that circulation forever. The understanding of 
The physics of circulation reduction is of major importance in the devel- 
opment of any phenomonological model and is discussed in greater detail in 
the following section. | 

2. Circulation Reduction 

All investigators using the DVM have recognized that greater atten- 
tion must be paid to the circulation reduction if the DVM is to produce 
results comparable to those obtained experimentally. In the early stages 
Of the evolution of the present model the need for a circulation reduction 
mechanism was quickly reaffirmed from the magnitudes of the lift and drag 
forces. 

Before undertaking a detailed discussion of possible circulation 
reduction mechanisms, it is advantageous to discuss the results obtained 
through the measurement of the vortex strengths. Schmidt and Tilmann [70], 
Through the use of an ingenious ultrasonic sound beam technique, measured 
the variations of the circulations of vortices without inserting a probe. 
Their measurements were made at various Reynolds numbers by changing both 
the velocity of flow and the diameter of the cylinder. Figure 15 shows T 
as a function of x and is representative of the results obtained by Schmidt 
and Tilmann. In their experiments the Strouhal number was about 0.21 and 


The mean velocity at the point of separation may be taken as U. - 1.5 U, 
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Fig. 15. Vortex strength versus distance [69, 701. 
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Fig. 16. Percent of circulation reduction versus distance E39, 20]: 
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Thus, the total circulation generated during the period T = I/f is approxi- 


mated by 
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In the observed formation region, | < x < 6, the vortices grow rapidly and 
at the same time lose some of their strength due to various circulation re- 
duction mechanisms. One can think of a net circulation reduction only for 
vortices which are in the stable and unstable regions. It must be noted 
that the ultrasonic technique used in Schmidt and Tilmann's experiments 
cannot yield accurate vortex strengths in the formation region because the 
path of the sound beam crosses over two oppositely signed vortices. In 
The stable region, however, only one vortex crosses the path of the sound 
beam at a given time. For the reasons cited above, part of The curve ex- 
tending from x = 10 to its experimentally reported maximum value of 40 is 
used to calculate the circulation reduction by dividing the circulation 
values by 10.7. The result is shown in Fig. 16. Also shown in Fig. 16 
are two reduction ratios computed similarly from Bloor and Gerrard's [69] 
reported vortex strengths at x values of 12 and 20. The circulation re- 
duction shown in Fig. 16 should be interpreted only as a reflection of 
the change of circulation with distance rather than as definitive values 
Since such values depend on the Reynolds number, three-dimensionality ef- 
fects, and the methods of their evaluation. 

One or more physical mechanisms are needed to explain and to model 
The observed circulation reduction. Circulation may be reduced only by 
Oppositely signed circulation. This hypothesis is somewhat mechanistic 
Dut is easily adaptable to numerica! modelling. There is no question that 


in reality the entire mechanism may be related to the action of viscosity. 
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Leaving aside momentarily the complex relationship between the action of 
viscosity and circulation reduction, attention ¡s concentrated on three 
sources of circulation exchange. It has been noted by the users of the 

DVM that the circulation may be reduced by the oppositely-signed circula- 
tion generated on the afterbody (countervorticity). In accordance with 
this assumption, the discrete vortices which come closer than a preassigned 
distance to the cylinder and the real vortices which cross the cylinder 
boundary are annihilated. Experiments by Fage and Johansen LIO] have shown 
that only about a 10 percent reduction is brought about by the entrainment 
of fluid bearing circulation of opposite sign from the flow aft of the 
separation points. Numerical experiments to be discussed in another sec- 
tion of this work, as well as those done by others, show that the circula- 
tion generated on the downstream face of the cylinder can account for no 
more than about 1/5-th of the tota! reduction in circulation. This fact 
calls for other mechanisms which account for the remainder of the circula- 
tion loss. The second mechanism postulated for the exchange of oppositely- 
Signed circulation is an integral part of the shedding mechanism. As pro- 
posed by Gerrard L63], when the growing vortex acquires a sufficiently 
large strength, it draws the feeding shear layer of the other growing vor- 
tex across the wake, thereby inhibiting further growth of the vortex. 

The tail of the cut sheet is then absorbed by the still attached vortex on 
the opposite side of the wake. This mechanism, quite distinct from the 
effect of the boundary layer on the rear of the cylinder, further reduces 
the strength of a given vortex. However, numerical experiments again show 
that the reduction due to this mechanism alone is far from sufficient to 
bring the circulation down to the observed or measured values. Consequent- 
ly, a third mechanism must be postulated with which not only the vortices 


in the formation region but also those further downstream experience 
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additional circulation reduction. The origins of the third reduction mech- 
anism may be partly in the ingestion or entrainment of fresh fluid into the 
Shear layers, partly in the instabilities in the spiralling vortex sheets, 
partly in slow viscous reduction, and partly through some other as yet 
unknown mechanism. Regardless of the true origins of the reduction mech- 
anism, two facts remain undisputed: the circulation of a given vortex does 
decay with time and the vorticity diffuses through viscous and Turbulent 
action; and, the DVM cannot predict even the most gross features of the 
real flows (for example, lift and drag coefficients and circulation distri- 
butions) without the use of a circulation reduction mechanism. 

lastly, it is necessary to briefly discuss the effect of diffusion 
of vorticity on the apparent reduction of circulation. Diffusion means 
That the vorticity is spread over a Bes: area. This in turn means that 
larger numbers of suitably distributed point vortices for each segment of 
the vortex sheet rather than a single elemental vortex represent the corre- 
ending distributed vorticity. Ordinarily, the point vortices (whether 
On a sheet or concentrated) are assumed to be placed at the center of 
gravity of the vorticity field.” This does not imply that the single vor- 
tex at the center of gravity will yield the same velocity at an arbitrary 
point on the cylinder or at the position of another point vortex as the 


distributed vorticity. As far as the cylinder ¡s concerned, the said 





The center of vorticity of M point vortices is computed as 
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effect becomes relatively small as the distance between the vortex and the 
cylinder increases. Consequently, the calculations are not noticeably 
affected by the diffusion of vorticity in the far field. Since diffusion 

is a matter of time, the vortices in the immediate vicinity of the cylinder 
do not undergo significant diffusion. The foregoing simply points out the 
difficulty of accounting for all direct or indirect decay mechanisms and 

the need for a suitable hypothesis which will account for the circulation 
decay and yet be compatible with the basic premises of the DVM. As is 
discussed in greater detail later, the present model accounts for the effect 
Of rear shear layer, circulation cancellation brought about by the shedding 
mechanism, and finally by the cross-street circulation exchange. The last 
process is accomplished simply by transferring a certain percentage circuia- 


tion between oppositely signed vortices at each time step. 


7.6 








IV. DETAILS OF THE NUMERICAL MODEL AND DISCUSSION OF RESULTS 


MEN SPECIFIC DETAILS OF THE MODEL : 
fe typical Sequence ot Calculations at an Arbitrary Time 

The most straightforward way to describe the particulars of the 
mode! is to examine the sequence of events that occur during a typical com- 
putational cycle. Figure I7 is a flowchart indicating the order in which 
The model executes the various phases of the analysis. Each of the events 
That occur during every computation cycle are described in the order of 
their occurrence. The asymmetry introduction and vortex shedding processes 
occur only under certain flow conditions and are addressed separately. 

At the beginning of each computation cycle the velocity distribu- 
tion on the forebody of the cylinder is calculated at one degree intervals 
from the most recent position of the forward stagnation point. The posi- 
tions of the upper and lower separation points are then determined through 
the use of Poh!hausen's method as described earlier. This allows the 
strength and positions of the two nascent vortices to be determined accord- 
ing to the no-slip condition. 

The next step is the convection of each of the vortices by a simple 
Eulerian method where the velocity of each vortex is computed at its pres- 
ent location. The velocities of the nascent vortices are computed at the 
point of their introduction prior to their actual placement into the flow. 
This insures that a nascent vortex is convected with the velocity prevail- 
ing at the time and place of its creation. Subsequently, the lift and drag 


coetficients are determined through the use of Eq. (15). 
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Vorticity Generation and Introduction 
a. Compute velocity distribution 
b. Determine separation points 
C. Assign nascent vortex strengths 
and positions 


Convection of Vorticity 
a. Compute vortex velocities 
D. Convect vortices 





Compute Lift and Drag Coefficients 


No hed 
3 
Vortex? 
Yes 
Cut Vortex Sheet 


| Rediscretize Vortex Sheets 
Circulation Reduct ion 


Display Results 


Fig. 17. Flowchart cf computation sequence. 
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The normal sequence of the calculations requires a test for the 
conditions to determine whether the vortex sheet should be cut and the 
attached vortex spiral be shed. Because of its importance and the need 
for special handling, it is discussed separately. | 

In discussing the problems associated with the use of discrete vor- 
Tices, it has been pointed out that the elemental vortices must be placed 
at the center of the sheet segment which they are supposed to represent. 
This led to the method of rediscretization. Thus, following the test for 
vortex shedding and if necessary, after the cutting of the vortex sheet, 
all attached and detached sheets are rediscretized. In doing so, each sheet 
is separately rediscretized (see Fig. 18). The discretization of the de- 
tached sheets is handled in accordance with the steps described in Chapter 
li, Section E. The only exception to the said process is In the treatment 
of the attached sheets near the nascent vortex. If has been repeatedly 
Pointed out that the conditions near the separation point must be disturbed 
as little as possible by rediscretization or any other flow of events. In 
accordance with this principle, the nascent vortex and its immediate 
neighbor are not allowed to change their strengths or positions. 

The next step in the sequence of calculations calls for circulation 
reduction and coalescence. It is simply hypothesized that every vortex in 
the flow field loses its strength in an amount proportional to its current 
strength. This hypothesis is based on heuristic reasoning, observations, 
and numerous calculations. Further justification of the hypothesis is made 
later. Here only the actual process incorporated into the computer code is 
described. 

During each computation cycle for t > 5, the strength of every vor- 


Tex, except the nascent vortices, is reduced according to the following 
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expression 
P (new value) = (l-p/100) T Cold value) (39) 


where p depends on the location KDE Reiz} as depicted in Fig. 19. For 
t< 5, p is varied linearly from zero to its final value of unity. 


Ip 
1.0 


Reiz} 
% 5 10 15 20 


Fig. 19. Variation of p with distance for t > 5. 


As discussed in connection with various circulation reduction mech- 
anisms, there are two other processes whereby the circulation ¡s reduced. 
The first is the reduction resulting from the proximity of vortices to the 
cylinder. Whenever a point vortex comes closer to the cylinder than a 
radial distance smaller than Ar = 0.04, that point vortex is assumed to be 
annihilated and removed from the flow field. The part of the sheet con- 
nected to the nascent vortex is excluded from this process. The annihila- 
Tion of circulation in this manner results only in a small reduction in the 
Total circulation and cannot account for a 40 to 50 percent reduction in 
The strength of the vortices in the stable region. The previous applica- 
tions of the DVM which considered this annihilation mechanism together with 
That resulting from the occasional proximity of oppositely-signed vortices 
(rediscretization not used) obtained only a mere 15 percent reduction in 
Eireulation [67, 681. This fact was the first indication of the need for 


waddi tional] circulation reduction mechanism. 
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The second of the two mechanisms concerns the entrainment of The 
tail of the detached sheet into the vortices across the wake. For the sake 
of discussion, assume that an attached vortex sheet is severed at a suit- 
able point. The details of how the cutting is accomplished is discussed 
later. The tail of the detached vortex sheet is drawn into the region be- 
tween two attached sheets because of induction. If left alone that tail 
eventually develops comp lex kinks and spirals and makes the rediscretiza- 
tion process unnecessarily complex. As proposed by Gerrard [63], the cir- 
culation in the tail gradually ends up in the attached sheet across the 
wake. In the model a similar procedure is employed. That portion of The 
tail which penetrates into the mixing region beyond a line tangent to the 
two attached vortex spirals (upstream of line AA' shown in Fig. 20) is re- 
moved from the flow field, and its circulation is uniformly subtracted 
from that of the sheet across the wake. Calculations which keep track of 
This cross-wake circulation reduction show that approximately five percent 
of the circulation of the detached sheet is transferred across the wake 
by this mechanism. It should also be noted that the position of the line 
severing the tail of the detached sheet is not of particular significance. 
One may have used other suitable and defensible procedures to accomplish 
the same purpose. 

2. Initiation of Asymmetry and the Vortex Shedding Processes 

The first of these two events occurs only once and the second every 
Time period a vortex is shed. As such, they constitute special events in 
The computational sequence and are discussed separately. 

a. Asymmetry Introduction 

In nature the inception of asymmetry and the shedding of vor- 
Tices are intrinsic characteristics of flow about bluff bodies. The causes 


of asymmetry and the subsequent vortex shedding are closely linked to the 
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instabilities in the shear layers. These instabilities may be brought about 
by the slight nonuniformities in the ambient flow, randomly distributed tur- 
bulence, asymmetry of the body, etc. It appears that The symmetric state 

of separated flow is not its most stable state. 

Whether the flow analysis is by the DVM, finite difference, 
finite element or whatever, asymmetry must be artificially introduced. 
Otherwise, the flow field remains forever perfectly symmetric and the vortex 
Shedding process never begins. Figure 2! illustrates the results of the 
present mode! when asymmetry is never forced to develop. 

Although the onset of asymmetry and the subsequent evolution 
of the wake have been observed and recorded by many researchers, little is 
understood about the controlling mechanism. The question of whether the 
asymmetry begins in the wake and is fed back to the boundary layer, or 
Vice versa, is unknown. Perhaps both occur simultaneously. Regardless, 

The process is not wel! enough understood to suggest any preferable way 

to artificially initiate its occurrence. Consequently, the method utilized 
in the present model is a result of heuristic reasoning based on wel | 
Known experimental observations and numerous numerical experiments with 
The model. Before describing the details ead results of the various af- 
tempts, it is worthwhile at this point to relate some of the general ob- 
servations, or lessons learned, from the investigation. Firstly, whatever 
method of asymmetry introduction is chosen, the steady state results of 
the model are identical; i.e., the model has no memory of the conditicns 
which initiated the shedding of the first vortex. Secondly, as stated 
earlier, asymmetry need be introduced only for the shedding of the first 
vortex. Following the shedding of the first vortex, the shedding process 
continues under its own impetus. Thirdly, every asymmetry mechanism in- 


volves a trade-off between the intensity and the duration of the perturbation 
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applied. Associated with the question of duration is the decision of when 
to initiate the mechanism. Fourthly, there is an optimum time during 
which the symmetric flow is most susceptible to the onset of asymmetry. 

If asymmetry is introduced during this period, the intensity and duration 
of its application can be minimized. Finally, of the unlimited number of 
mechanisms that could be devised, it is best to choose the simplest method 
that produces the desired asymmetry. 

The interactive nature of the vorticity generation mechanism 
and the vorticity distribution in the wake suggest that asymmetry might be 
Introduced either into the generation mechanism or into the wake itself. 
Methods applied to the vorticity generation method involve, for example, 
artificially increasing the computed vorticity flux at one separation point 
(say by 10 percent) and simultaneously decreasing the flux at the other 
Separation point. However, sample calculations with this procedure show 
that the method is not effective. Furthermore, the large asymmetry causes 
the computed separation angles to undergo unrealistic excursions from their 
mean positions. 

Introduction of asymmetry into the wake can be accompiished by 
either artificially altering the strength and/or position of the point vor- 
tices on the existing sheets. Changing the strengths of the vortices in 
one sheet relative to the other results in a sluggish response of the wake 
Similar to that described above. The general ineffectiveness of altering 
vortex strengths to produce the downstream convection of one sheet relative 
To the other suggests the more direct approach used in the present model, 
The rationale [sto initiate the downstream convection of one sheet by 
actually displacing it a smal! amount each computation cycle over a short 
Period of time so that when the mechanism is discontinued the sheet con- 


tinues to move without external influence. 
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Numerous numerical experiments have shown that a gradual appli- 
cation and buildup of the displacement increment serves to provide a smooth 
transition from symmetric to asymmetric conditions. The initiation time 
and duration of application are based primarily on two considerations. 
Firstly, this time interval should correspond to the period during which 
the flow is most susceptible to asymmetry. Secondly, the initiation of the 
mechanism should begin after the initial peak in the drag curve. Since the 
cause of this drag overshoot is the AT growth of the vortices, the 
asymmetry mechanism should not be allowed to influence this symmetric de- 
velopment. 

In the following the specifics of the asymmetry introduction 
method are presented. During each computation cycle within the interval 
9.0 < t < 9.0 the vortices of the upper attached sheet are displaced down- 


Stream an amount given by (see Fig. 22) 


Ax = 0.08 {1.0 - cos[m(t-5)/4]} (40) 
016 [4 
0.08 
SES 5 6 7 8 9 


Fig. 22. Variation of the asymmetry perturbation Ax with time. 
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Several points are worthy of discussion. Firstly, the displacement given 
by Eq. (40) is not unique and any other form of displacement function could 
have accomplished the same purpose. The particular choice, depicted by 
Fig. 22, introduces the disturbances gradually and continuously without 
shocking the flow field. Lastly, the total displacement applied to the 
emcer during the time interval 5 < f < 9 is only 32 percent of the radius 
of the cylinder. 

b. Vortex Shedding Process 

Vortex shedding is understood as the mechanism whereby the feed- 

ing layer is cut. The specific details of the conditions which dictate the 
separation of a spiralling sheet from its feeding layer are not known and 
have yet to be determined experimentally. The numerical model is not in 
need of such a cut were it not due to the limitations of computation time. 
If the flow field is allowed to evolve asymmetrically, the flow pattern 
shown in Fig. 23 results. This figure shows not only the evolution of the 
doubl y-connected vortex sheets but also the drawing-in of the tail of the 
shedding vortex into the mixing region. The continuation of the calculations 
in this manner are prohibitive even with the largest computers. However, 
it is instructive in two ways to examine the evolution of the wake in this 
Manner. Firstly, the uncut sheet begins to change the sign of its curva- 
Ture as the lift force goes through its extremal (see point A of Fig. 23). 
The change of curvature is in fact the beginning of the drawing-in process 
of the shear layer across the wake. Subsequently, a strong inflection and 
Kink develop drawing the sheet closer into the formation region. Simul- 
taneously, the lift force increases towards its next extremal. The gradual! 
nature of this process does not permit one to identify a specific moment 
at which the shedding should commence. Thus, the examination of other flow 


features exhibiting stronger reaction to the precise moment of shedding are 
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necessary. Numerous calculations with symmetric and asymmetric evolution 
of the wake and with and without the rediscretization of the shear layers 
have shown that the sheet fo be cut bears its minimum circulation per unit 
length at a time immediately after dT/dt goes through a minimum (see Fig. 
24). In Eher words, the connecting sheet is weakest and, therefore, 
most susceptible to breaking when ifs circulation is minimum. IT is on 
the basis of this fact that the sheet is cut immediately after the rate 
of circulation reaches its minimum. At the time of cut, the lift force 
is about 50-percent of its preceding extremal (see Fig. 25). 

The actual cutting process is accomplished by removing a single 
point vortex, for only one time step, from The sheet at a distance s = 0.4 
(Fig. 26). Subsequent calculations separately rediscretize the detached 
Sheet and the part still connected to the nascent vortex. As The flow 
evolves, the attached sheet continues To grow rapidly to form a new spiral 
while the detached sheet moves further downstream. AT each Time step the 
computer code tests the condition as to whether The new connected sheet 
should be cut. At The time The said condition is satisfied, the vortex 
sheet is cut as described above. It is clear that the shedding of the vor- 
Tices and the cutting process alternate with a time period equal to half 
The Strouhal period. When a vortex sheet is cut, say at the lower side of 
the cylinder, the sheet which was previously detached from the top part of 
the cylinder is coalesced into a single vortex at its center of gravity. 
At the time of its coalescence, the core of the vortex spiral has moved a 
distance at least seven radii downstream. The calculations have conclu- 
Sively shown that the process of coalescence at this distance from the 
cylinder does not alter the continuity of any variable (e.g., lift, drag, 


Circulation, etc.). 
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Cutting of a vortex sheet. 
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DISCUSS ON OF RESULTS AND SENSITIVITY ANALYSIS 
in the following, the results obtained with the numerical model (re- 
ferred to hereafter as the standard run) are presented. Following this, 
the variation of some of the important kinematic and dynamic parameters 
to the variation in others is discussed. The purpose of these discussions 
is to partly establish the stability of the results, partly to relax the 
range of variability of the nondisposable parameters, and partly to estab- 
lish the cause and effect relationships. AT this time, no distinction 
is offered concerning the choice of parameters in the model and the calcu- 
lated results. Furthermore, no distinction is made between prediction and 
inspired hindsight. 
I. Results Obtained with the Standard Run 

The standard run is computed using At = 0.125, p = 1.0 (see Fig. 19) 
and with the asymmetry introduction and the shedding mechanisms previously 
described. 

The evolution of the wake over one arbitrarily chosen cycle is shown 
in Figs. 27 and 28. Figure 27 depicts the near wake region and Fig. 28 
shows the far wake vortex street. 

The variations of the lift and drag coefficients with t are shown 
in Figs. 29 and 30. The drag coefficient reaches a maximum at + = 4.2 and 
then rapidly decreases to its ultimate mean value of about 1.2. The fluc- 
tuations of Co have a frequency twice the vortex shedding frequency and an 
amplitude of about 0.06. The overshoot of the drag coefficient is a direct 
consequence of the rapid accumulation of vorticity in the symmetrically 
growing vortex spirals. Following the shedding of the first vortex, Ch 
drops sharply and the lift coefficient increases, ultimately reaching a mean 
amplitude of about 0.65. The slight variations in the oscillations of C 


D 


and CL come from various sources. The vortices absorbed during the cross- 
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Fig. 27. Evolution of flow in the near wake. 


95 








Rø i 
JOD 


QRS 


Fig. 27 (con't). Evolution of flow in the near wake. 
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Fig. 27 (con't). Evolution of flow in the near wake. 
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wake cancellation mechanism may be one more or one less during two different 
cycles. The increase of the otherwise finite length of the wake does not 
permit the vortex sheets to follow exactly identical paths from one cycle 
to another. Regardless, the said perturbations in Co and Ch are quite 
small and well within the accuracy of the numerical model. 

The amplitude of the lift coefficient is certainly within the 
range of lift coefficients reported in the literature. However, if is not 
possible to find a truly two-dimensional laboratory experiment to compare 
The calculated results with greater definitiveness. More is said later 


about the dependence of the amplitude of C, on other parameters, in par- 


L 
Ticular on the circulation reduction mechanism. 

Figure 3| shows a comparison of the calculated and measured drag 
coefficients. The results are strikingly similar and may be taken as a 
first indication of the workings of the model. 

The oscillations of the separation points in terms of the angles 
measured from 8 = O (not necessarily the instantaneous position of the stag- 
nation point) are shown in Fig. 32. The separation angle for a given sepa- 
ration point oscillates nearly sinusoidally with an amplitude of A8 = 5 
degrees about a mean value of 8. = 17 degrees. The frequency of oscilla- 
mon of an is identical to that of the shedding of vortices. 

The stagnation point does not remain at 9 = O even in steady flow 
about a stationary cylinder. Figure 33 shows the oscillation of the stagna- 
tion point. Similar oscillations have been reported by Dwyer and McCroskey 
[71]. The separation angles relative to the instantaneous position of the 
Stagnation point are shown in Fig. 34. A comparison of Figs. 32 and 33 
Show that the oscillations of the separation and stagnation points are 180 
degrees out of phase. In other words, if the stagnation pofnt moves up, 


say to 9 = 2, the separation point at the upper half of the cylinder moves 
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backwards from a mean position In = // to a. = 74 degrees. This fact is 
also in conformity with previously reported results [71]. 

The rate of change of circulation at both the upper and lower sepa- 
ration points is shown in Fig. 35. Because the generation and shedding of 
Circulation are the essence of the entire bluff body problem, dI'/dt has a 
special significance among al! other results predicted by the numerical 
model. A cursory examination of this figure shows that there is a perfect 
Similarity of the rates of circulation of the upper and lower nascent vor- 
tices. Secondly, the variations in dP/dt from one cycle to another as t 
increases from 40 to 200 are very minor, indicating that the vortex street 
may be regarded as Infinitely long for all intents and purposes for t 
larger than about 40. 

The mean value of |dT/dt| is about 1.2. Noting that the mean base 


pressure C is related to |dT/dt| by [72] 


pb 


Cc, = 1.0 ~ 2 [df /at| (41) 
pb 


The present calculations yielded Sote = -|.55 which compares favorably with 
that reported by Roshko LII_] at high subcritical Reynolds numbers [Cb 
(experimental) = -1.36 for 20,000 < Re < 100,000]. 

The amplitude of oscillation of dl'/dt is about 0.18 which indicates 
That the velocity at the separation point varies from |.4! to 1.65. 

All of the parameters cited above,namely lift and drag coefficients, 
Separation and stagnation angles, and dI'/dt oscillate with a frequency 
identical to the vortex shedding frequency f. The Strouhal number, St = 
f D/U is found to be St = 0.205. The experimental values of the Strouhal 
number are reported to lie between 0.195 and 0.210 for the Reynolds number 


range from 2,000 to 400,000 (see Fig. 36). Thus, the predictions of the 


model are in conformity with the measurements in the range of Reynolds 
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Fig. 36. Strouhal number versus Reynolds number [73]. 





numbers where the DVM is applicable. As is noted later, the Strouhal number 
is the least sensitive of all the parameters to the variations of the non- 
disposable parameters. 

One of the often discussed characteristics of The vortex street, 
particularly in the stable region, is the relative spacing of the vortices. 
Yon Karman's classic stability analysis gave WEN = 0,281 (see Fig. 37). 
Since von Karman, the said ratio has been determined experimentally by a 
great number of investigators and it has been found That hR, varies from 
Braut 0.19 To 0.5. 

The lateral and transverse vortex spacings obtained with the 
Standard run are shown in Fig. 37. The asymptotic values of N and h, are 
well within the experimentally determined values [74]. In the stable region 
of the street, the lateral spacing decreases somewhat and the transverse 
Spacing gradually increases from about 0.63 to its asymptotic value of about 
E05. The ratio AA, in the asymptotic range is seen to be 0.23, a value 
which compares well with those reported experimentally [74]. Evidently, 
ir is not meaningful to think of vortex spacings within the formation 
region. 

All discrete vortex models of separated flow about bluff bodies 
indicate the need for greater loss of circulation in the wake in order to 
bring a closer correspondence between the calculated and measured separa- 
tion points, lift and drag coefficients, dl /dt, etc. For example, Deffen- 
baugh and Marshall [47] obtained a mean separation angle of 67 degrees and 
a mean drag coefficientof 0.9 even though the smaller separation angle 
should have yielded larger drag than that found experimentally due to the 
increased wake size. In the present model all conceivable circulation 
reduction mechanisms have been taken into consideration. Their variations 


with x or t may not be unique or in conformity with those obtained 
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experimentally but their characteristics follow all vortex shedding mech- 
anisms so far advanced. Figure 38 shows the circulation retained by a 
vortex sheet from its inception to its subsequent coalescence and convec- 
tion into the far wake. Also shown in this figure are the two experimen- 
tal points reported by Bloor and Gerrard [69] and the mean curve drawn 
through the data reported by Schmidt and Tilmann [70] (see Fig. 15). The 
comparison of all the results is better than expected in view of the fact 
that the experimental determination of circulation through direct or 
indirect methods is an extremely difficult and approximate task. In par- 
ticular, the ultrasonic sound beam technique used by Schmidt and Tilmann 
yields the vortex strengths only approximately in the stable and unstable 
region of the vortex street. There are no measurements of the vortex 
strengths in the formation region. The foregoing does not imply that the 
particular variation assigned to the circulation loss mechanism (see Fig. 
19) is unique. In fact, the only major nondisposable parameter in the 
entire numerical model is the form assigned to the circulation loss. The 
need for circulation loss has been amply demonstrated. The form of the 
particular er Rit can only be adjusted with hindsight by comparing 
all of the predictions of the model with those obtained experimentally. 
Because of the fact that the differences in measured and calculated circu- 
lations are further downstream in the wake, the use of the experimental ly 
suggested p variation does not in any way affect tift, drag, separation 
angles, dP/dt, and the Strouhal number. This wil! become more evident in 
The course of the discussion of the sensitivity of the model. Calculations 
not reported herein with more exotic forms of p have ascertained the 
validity of the insignificance of variation of p beyond the stable region. 

2. Sensitivity Analysis 

This section deals with the effect of the variation of one parameter 


abouts its value used in the standard run on the remaining major parameters. 
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All numerical analyses must use a finite mesh size or time interval. 
The magnitude of these quantities might affect the results in various ways. 
In the DVM a large At might bring in curvature effects particularly in flow 
regions where the radius of curvature of the spiral is small. This effect 
is brought about by the convection of the vortices along straight lines co- 
incident with their instantaneous velocities. Thus At should be kept reason- 
ably small to allow the vortices to follow their natural paths. A very 
small At, on the other hand, requires prohibitively large computation times 
and defeats the purpose of the numerical ee 

The standard run was made with At = 0.125. Additional runs were 
made varying At from 0.| to 0.2. The resulting drag and lift coefficients 
are presented in Figs. 39 and 40. The differences in the two coefficients 
for all values are negligibly smal! and certainly within the range of 
accuracy expected from numerical calculations. Other parameters such as 
separation and stagnation angles, circulation loss, Strouhal number, etc. 
remain essentially identical to their corresponding values obtained with 
Bp = 0.125. 

The role played by the loss of circulation in both laboratory and 
numerical experiments has already been discussed. Here attention is devoted 
to the selection of the particular value of p and the effect of its varia- 
tion on the remainder of the predictions. Several hundred hours of calcu- 
lations have shown that one can arrive at an approximate value of p by 
comparing the predicted and measured values of the vortex strengths and the 
lift and drag coefficients. None of these parameters, however, exhibit the 
Same degree of variation with p at corresponding times and downstream dis- 
Tances. For example, the comparison of the vortex strengths in the stable 


region, where most of the measurements are made, is not very meaningful 
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Fig. 39. Drag coefficient versus time for At=0.1, Oras ad Oa: 
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since the effect of such vortices on the lift and drag coefficients is very 
small. Thus, it was decided to choose the p value such that the terminal 
value of the drag coefficient conforms to that measured in the subcritical 
region. No other experimental data was used in The analysis. Numerical 
experiments with various values of p have shown that p = |.0 yields Co = 
1.2. Having determined a value of p, it is necessary to examine its varia- 
tion on all of the predictions of the model. 

The drag coefficients obtained with p = 0, p = 0.5, p = 1, and 
p = 2 are shown in Fig. 41. The insert in Fig. 41 shows the variation of 
Co with p in the asymptotic range of Ch. 


p about p = | caused only about 10 percent change in C 


Evidently a 100 percent change in 


9° Consequently, the 


Tuning of the model through the use of Ch = AS MON As Critical as it 
may have appeared at first sight. Figure 4l also shows that a circulation 


loss mechanism is necessary (see curve for p = 0) if the numerical model 
is to yield results in conformity with the experiments. As noted earlier, 
the need for circulation loss is also evidenced by the vortex-strength 
measurements (see Fig. 38). Thus, the important question in the evolution 
of the model is not whether there should be a circulation loss but rather 
what should the mechanism and magnitude of this loss be. 

The variation of the lift coefficient with p is shown in Fig. 42. 
It is immediately apparent that the lift coefficient is far more sensitive 


Kè p than C Regardless of its sensitivity, however, C, remains well 


p' L 
within the experimentally reported values. One could fine-tune the value 


of p by matching the measured and calculated values of C Unfortunately, 


L” 


fms 15 not possible since there are no reliable €, values, particularly in 


E 
the region of high subcritical Reynolds numbers. Consequently, no attempt 


is made to slightly increase or decrease p from p = | to obtain a specific 
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Fig. 41. Drag coefficient versus time for p=0.0, 0.5, |.0 and 2.0. 
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CL value. The importance of Fig. 42, rather, lies in the fact that the 


numerical model clearly shows the sensitivity of C, to circulation loss and 


L 
offers an indirect explanation of the difficulties encountered in the 
laboratory experiments. The findings of the mode! can also be taken as a 
means of reducing the transverse force by artificially increasing the 
circulation dissipation in vortices (e.g., use of fairings, wire screens, 
eric.) . 

The Strouhal number displays a relatively small variation with p 
(see Fig. 43). This explains the good correspondence obtained between 


measured and calculated values of St even when the circulation reduction 


does not exceed |5 percent (Clements L67], Sarpkaya [68]). 
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Fig. 43. Variation of the Strouhal number with p. 


It is of Interest to point out the relationship between the strength 
and motion of vortices in the near wake and the smal! variations in the 
Strouhal number. When the detached vortex is relatively stronger (p < I), 
then it lingers longer near the cylinder by virtue of the velocities in- 
duced on it by the equally-strong vortex across the wake. This in turn in- 
creases the period of vortex shedding or decreases the Strouhal number. 


For weaker vortices (p > |) the inverse is true. 
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In the following, the sensitivity of the lift and drag coefficients 
and the Strouhal number to the variations in the amplitude of p has been 
discussed. In doing so, the shape of the dissipation function has not 
been changed. From a scientific point of view, the next to consider is 
The effect of the variation of the dissipation function on the major predic- 
tions of the model. For this purpose an extensive study was undertaken To 
acquire a deeper understanding of the dissipation particularly in the 
near wake region of the cylinder. In doing so, not only the terminal 
value of the drag coefficient but also other parameters such as fhe lift 
coefficient, separation points, and the amount of circulation retained by 
the concentrated vortices were carefully observed. Here only the high- 
lights of these observations can be described. The use of a dissipation 
function such as p = o/x (for x > |) where o is a constant shifts the re- 
gion of larger dissipation closer to the rear of the cylinder. For example, 
The results obtained with o = 2.0 have shown that the lift and drag coeffi- 
Cients show no more than a 10 percent variation from those obtained with 
the standard model. The separation points remained at almost exactly the 
Same average values. The percent circulation retained by the concentrated 
vortices increased for obvious reasons. In fact, the relative vortex 
strengths at comparable distances were somewhat larger than those measured 
by Gerrard [24]. The Strouhal number decreased from 0.205 to 0.193, in- 
dicating that not only the vortices ¡in the formation region but also those 
further downstream determine the transport velocity of the vortices and 
hence the Strouhal number. 

The experience gained from this exploratory investigation has shown 
That there should be a fairly large and uniform dissipation in the near 


wake, say within the formation region, and then decrease gradually with 


ars 





distance. However, without careful measurements of the vortex strengths 
as a function of time or distance it may never be possible to determine 
the precise shape of the dissipation function. 

To further determine the sensitivity of the model to dissipation and 
The relationship between some of the most important parameters (e.g., 
separation angles, vortex strengths, and Strouhal number), an additional 
Series of numerical experiments have been performed. Specifically, The 
Separation angle was determined not in the manner described in the stan- 
dard model but rather by suitably extrapolating the velocity profiles 
beyond the point of maximum velocity. The use of such an extrapolation 
alleviated the influence of the nascent-vortex proximity effect on the 
Separation point and resulted in larger separation angles (8. = 8i degrees). 
When the shape of the dissipation function was kept exactly as in the 
standard run, the drag coefficient reached an ultimate value of about 1.5 
and the lift coefficient a value of about 0.8. The Strouhal number de- 
creased to about 0.19. [tT was thus apparent that the value of p in the 
Standard dissipation should be somewhat scada Say tò prè 1.5.9 Suen 
an increase in p would decrease both Co and SL and increase the Strouhal 
number without affecting the separation and stagnation points. This 
numerical experiment suggests that it is possible to match the observed 
and culata separation points and resistance to fluid motion with one 
nondisposable parameter. No attempt has been made to find a precise value 
or shape for the dissipation function because of the primary concern with 
The interrelationships between various parameters. The fact has been 
demonstrated that the physical mechanism governing the flow is not affected 
by the actual value of p or C. and C.. It is evident from the foregoing 


D E 


that the numerical modelling of a turbulent wake does not require a 
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dissipation based on turbulent diffusion and that diffusion could be 
quantified only by additional measurements which will complement the 
numerical experiments. As is often the case In many aspects of fluid 
mechanics, the question seems to have been reduced once again to the 


understanding of the physics of turbulence. 


12 | 





V. SUMMARY OF THE RESULTS 


The numerical model based on the rediscretization of vortex sheets and 
The mutual interaction between the boundary layers and the wake yielded 
numerous results which may be summarized as follows: 
I. A discrete vortex model of flow past a stationary cylinder may be for- 
mulated with only one nondisposable parameter. Experiments and the numeri- 
cal results show the necessity of circulation reduction but not necessarily 


The exact mechanism by which it is accomplished. 


2. Based on a postulated circulation reduction mechanism, all predictions 
of the model are in conformity with those obtained experimentally. The 
fact that all and not only a few of the parameters agree with the experi- 
mental measurements adds further credibility to the assumptions made and 
to the power of prediction of the model. 

5. A concerted effort is made to establish a balance between simplicity 
and accuracy. To this end, simple linear relationships have been preferred 
in lieu of more complex ones. 

4. One of the most important aspects of the numerical model is its flexi- 
bility. Practically all parameters have a wide stable range and allow 

The user to perform numerical experiments to seek a physical understanding 
of the intricate relationships between two or more parameters. 

5. The demonstration of the applicability of the discrete vortex mode! 

To the analysis of flow about a bluff body with mobile separation points 
Shows that similar models can be constructed with ease for bodies with 
fixed separation points. This allows one to study such problems as vortex 


shedding from ship bilge keels, stall over aircraft, fuel slushing in large 
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containers, etc., in naval hydrodynamics and aeronautics. 

6. The model, as it stands, can be used in the analysis of a great many 
flow-structure interaction problems. Of particular interest is the analysis 
of a circular cylinder undergoing forced or self-excited transverse oscil- 
lations. Clearly, the numerical experiments with oscillating cylinders 

can not only shed additional light on the mechanisms which bring about such 


oscillations but also increase the untversality of the model. 
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USER AP EA I ONEORTEEEZDISEREHEWIORHET MODEL 
TO EN EXCITED TRANSVERSE OSCILLATIONS 
A. INTRODUCTION : 

Numerous experiments have shown that (see Parkinson [75] for an ex- 
tensive review) when the natural frequency of a bluff body is close to 
the vortex-shedding frequency, The body may undergo self-excited trans- 
verse oscillations. The primary consequences of these oscillations are 
as follows: (i) The vortex shedding frequency locks on to the natural 
frequency of the body. (ii) The amplitudes of oscillation and the 
Transverse force increase rapidly up to the point of perfect synchroni- 
zation where the maximum amplitude is attained. Subsequently, the oscil- 
lations decrease either rapidly or abruptly. The abrupt drop in amplitude 
iS a consequence of hysteresis. The reasons leading to the hysteresis 
are not yet known. (iii) The transverse force leads the displacement 
by a phase angle which increases gradually at first to about 50 degrees 
and then rapidly to about 130 degrees. (iv) At the end of the synchroni- 
zation region, the vortex shedding frequency jumps to that governed by 
the Strouhal relationship. 

The studies concerned with the understanding and prediction of the 
synchronization phenomenon may be classified in three categories. The 
first and most extensive category has been the experimental observations 
and measurements with forced or self-excited oscillations of circular 
cylinders, cables, and a smal! group of other bluff bodies. The second 
category of investigations dealt primarily with the analysis of the 


Structural aspects of the equations governing the vibrations by making 
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ad hoc assumptions about driving fluid forces. Finally, the third cate- 
gory concentrated on the understanding of the flow field and the nature 
of the fluid forces. It has not yet been possible to understand the 
coupled structural and fluid mechanical mechanisms to explain the reasons 
leading to the self-excited oscillations. 

Mathematical models based on The assumption of a nonlinear oscillator 
have been developed. Among those, the most noteworthy is that proposed 
by Hartlen and Currie [76]. Their model predicts some of the observed 
features of the vortex-induced oscillations. However, the model requires 
a number of nondisposable parameters and gives no clues about the cause 
and effect relationships. IT has been realized, as suggested by 
Parkinson [75] among others, that if any understanding of the phenomenon 
is to be achieved, theoretical developments and numerical experiments 
must be compared with laboratory experiments to guide and complement each 
other. If is with this suggestion in mind that the discrete vortex model 
has been applied to the prediction of the characteristics of ee: 
oscillationsofa circular cylinder. In the following, first the 
description of the mathematical ar e of the problem, then the 
specific details of the numerical procedure, and finally the discussion 


of the results and a summary are presented. 


8. FORMULATION OF THE PROBLEM 
The equation of motion for an elastically-mounted and linearly-damped 
Cylinder of unit length (see Fig. 44) may be written in dimensional form 


as 


my + 2m W. cy +m 0 or Ue C, (42) 
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where m represents the mass of the cylinder; y, the displacement; Ws the 
natural circular frequency; Z, the damping ratio; Pes the fluid density; 
e, the radius of the cylinder; and Ch, the instantaneous value of the 
lift coefficient. In Eq. (42), y and Y represent respectively the first 


and the second derivatives of y with respect to Time. 





; 2 
Zmw CY mw. y 


Fig. 44. Elastically-mounted, linearly-damped cylinder. 


Defining fø as the vortex shedding frequency for the stationary 
O 
cylinder, Je as the corresponding Strouhal number, and 
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Eq. (42) may be written in normalized form as 
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ese a net: a n = 27 S= ac (44) 
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Because of the need to distinguish the stationary and oscillating vortex 
shedding frequencies, and hence the Strouhal numbers, the previously and 
common |y used notation of St for the Strouhal number is modified as 
denoted above. 

The evaluation of Eq. (44) requires the instantaneous value of The 
lift coefficient. In the present model it is evaluated through the use 
of the complex velocity potential modified to take into account the motion 


of the cylinder. Thus, one has 








> N 
[=z 
Der E22 )# ——H Im Y A eZ 2 = J} (45) 
n= | n oO 


in which ZO denotes the instantaneous position of the center of the cylin- 
der and u its velocity (see Fig. 45). The use of Eq. (45) together with 


the generalized Blasius theorem yields 


N 


D - my (46) 
n O 


n= 


8 pil, 





Fig. 45. Coordinate system for oscillating cylinder. 
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A comparison of Eqs. (46) and (15) shows that Z a is replaced by the rela- 
tive position 2 72, in the latter equation and two new terms are added 
which account for the forces arising from the velocity and acceleration 
of the cylinder. 

The calculations for the boundary layer and The separation points re- 
quire the evaluation of the tangential velocities relative to the cylinder. 
At an arbitrary time + and angle 9 (measured as shown in Fig. 46), the 


tangential velocity is given by 


-i (0-7/2); 








= Re {Latz_ +e 9) -z le (47) 
X 
Fig. 46. Relative tangential velocity on boundary of 
moving cylinder. 
As before, the velocity at any point in the flow field is given by 
|-ż : 
: dw “O i | 
=UPTIV = a— = m se — Da e (48) 
Te ee "om Lot n Z-Z (z-z_)- — 
O n= | o — 
2-2. 
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C. NUMERICAL PROCEDURE 

lt is apparent from the formulation of the problem that the solution 
of Eq. (44) requires two parameters (Z and a) characterizing the mechani- 
cal system and the instantaneous fluid force imparted to the cylinder. 
The parameter a, known as The mass parameter, may be readily calculated 
for a given cylinder of mass m, radius c, and fluid density O [see Eq. 
(43) ]. The damping ratio, Z, is a measure of the relative energy loss 
within the SA due to molecular motion and does not include the so 
called fluid damping. All fluid forces acting on the cylinder are included 


in the proper specification of C, which represents the integrated effect 


E 
of the surface-pressure loading by the fluid on the vibrating cylinder. 

The discrete vortex model developed for the stationary cylinder is 
used together with Eqs. (46) and (48) to determine the fluid forces acting 
on the cylinder. Equation (44) is integrated through the use of a fourth 
order Runge-Kutta method. Otherwise, no additional modifications are 
introduced into the procedures described in the previous chapter. In fact, 
the only specification made was the time at which the cylinder was al lowed 
to respond to the fluid forces. Such a decision was necessary to avoid 
the coupling of two transient states resulting from the impulsive start 
of the flow and the transient response of the cylinder. This would have 
made the calculations unnecessarily complex. Instead, it was decided 
first to establish a flow with a sufficient number of vortices in the 
wake (e.g., by regarding the cylinder as stationary for t < 10) and there- 
by studying the response of the cylinder through a more or less steady 
uniform flow. The approximations signified by "more or less" are due to 
The fact in all numerical calculations the vortex street is finite. In the 
present calculations, the cylinder is allowed to respond after t > 10, 


which was sufficiently large for all intents and purposes. 
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Mathematically speaking, there are no particular restrictions on the 
selection of the numerical values of Z and a. other than the fact that 
they should be physically realizable. From the point of view of actual 
calculation, however, one must place some restrictions on £ and a. TO 
limit the computation time. The smaller the value of Z, the larger the 
amp!itude of oscillation. The larger the mass of The cylinder or the 
smaller the parameter a. (keeping everything else constant), the larger 
is the transient Time during which The cylinder oscillations asymptoti- 
cally approach their final values. Consequently, a very small value of 


¢ coupled with a relatively small value of a, would require many cycles 


Ò 
of oscillation before the cylinder reaches its final amplitude. This 

makes the computation prohibitively expensive and time consuming. The 
foregoing arguments constituted one of the primary conditions in the selec- 
tion of Z and ase The second and probably more important reason from a 
scientific point of view was the focusing of attention on the understand- 
ing of the physics of the phenomenon rather than on the solution of 
specific problems. With the foregoing considerations in mind, the Z and 


a. values shown in Table II were used in the calculations. The role of 


the new parameter Se = c/a, appearing in the Table is discussed later. 


TABEEZUN 
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Be DISCUSSION OF RESULTS 
e ase aUe pen 

The major part of the discussion of The results deals with the 
Phenomenological predictions of the model and their comparison with those 
observed or measured experimentally. For this purpose only one set of 
damping ratio and mass parameter is considered. The variations of The 
lift coefficient, relative amplitude, phase angle between the lift force 
and the cylinder response, and finally the ratio of the cyl inder-osci | la- 
tion frequency to natural frequency is presented as a function of U 


defined by 


Oe (49) 





In dimensional terms, Un = U/w D . The reasons for the choice of this 
particular parameter ae twofold. Firstly, it is commonly referred To 
in the previous works [75]. Secondly, it is a measure of the velocity 
of flow past a cylinder of given diameter and natural frequency. 

Figure (47) shows the results obtained by Feng [77] with a lightly 
damped circular cylinder. As noted earlier, when the vortex shedding 
frequency, Te , for a stationary cylinder approaches the natural frequency 
EU = 0.8 in Fig. (47) ], the vortex shedding frequency becomes nearly 
identical to the natural frequency of the cylinder. Furthermore, the 
frequency of oscillation of the cylinder, Wo, nearly coincides with ¡Ts 
natural frequency. As U increases towards |.O, the amplitude of both 
the oscillation and the transverse force increases rapidly while w/w and 
w/a remain nearly constant at a value slightly under unity. As U in- 
creases further, both amplitudes drop either gradually or abruptly as 


illustrated in Fig. (47). The abrupt drop is a consequence of hysteresis 
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Experimental results for self-excited transverse oscilla- 
tions of a cylinder [77]. 
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noted earlier. For U > |.| the frequency of vortex shedding, W,» jumps 


to its stationary-cylinder value, i.e., w, 3w. However, the cylinder 
continues to oscillate at the natural A ae the cylinder with very 
small amplitudes. One can, therefore, identify two major regions in the 
relationship between My Us and We For U smaller than about 0.8, w, = 


W and w = Pine ne rande OOc U -cs TNN w L and» zu 2u- 
‘a Oa We 9 a > ES: y VE V € n 


< 195 oO, = My and W = W, e 

The phase angle between the exciting force and the cylinder re- 
Sponse undergoes dramatic changes in the region where 0.8 e U < PEP 
The reasons leading to the amplification of the oscillations of the 
cylinder and the lift force, and the rapid increase of the phase angle 
from about 20 degrees to 100 degrees have not been understood because of 
The difficulty of measuring or observing The detailed kinematics of the 
flow field about the cylinder. Many experiments have essentially confirmed 
the general features of The results obtained by Feng E77]. “Most of The 
attention has been devoted to The understanding of the locking-on of the 
vortex shedding frequency to the natural frequency of the cylinder in the 


range of U values 0.8 <U_<l.l. In conformity with the previous 


e 

literature, this region is variously referred to as the synchronization 

region, locking-on region, capture region, or simply as the resonant region. 
The numerical results obtained with Z = 0.02 and a ae 0.0125 are 

shown in Fig. (48), again as a function of U. It is immediately apparent 

that the results shown in this figure follow essentially the same trends 

as those in Fig. (47). Also shown in Fig. (48) are the corresponding 

values of w and ratios w / wr or U ZUÈ where ws and Uf are the values of 


W and U when Ny reaches its maximum value (point A). This point is 


referred to as the point of perfect synchronization. 
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Fig. 48. Numerical results for self-excited transverse 
oscillations with 2=0.02 and a 70.0125. 
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2. Mechanism of Synchronization 


Previous investigators 75] have uncovered a number of relatively 
Unrelated facts regarding the kinematic and dynamic behavior of flow 
about a +transversely-osci | lating cylinder. The most notable of the ob- 


Servations and measurements may be summarized as follows: 


mi Oscillations increase the spanwise correlation; 

eit) Miewsimengiaot ane YORI ices Increases; 

(iii) The mean base pressure decreases and the mean drag coefficient 
increases; 

iV) Transverse force is amplified relative to its stationary-cylinder 
value; 

(v) The phase angle between lift and cylinder motion increases with 
Un 

(vi) The velocity range over which synchronization occurs increases 


with oscillation amplitude; 
(vii) Separation point excursion increases with oscillation; 
(viii) Lateral spacing of the vortices decreases with increasing ampli- 
tude and is unaffected by changes in frequency of oscillation; 
CTx) Wake-formation length decreases systematically with increasing 
amplitude of vibration for U < ee and increases for U > UT; 
(x) The wake contracts for U.< Ue and expands for ve ur. 

The facts summarized above have been distilled from a number of 
investigations, not all of which have been conducted at identical Reynolds 
numbers or under similar experimental conditions. Consequently, there 
are certain features of the flow which have been observed to increase in 
one investigation and decrease in another. Whenever such conflicts arose, 


heavier weight was subjectively attached to the observation more frequently 
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reported. Perhaps the most ¡nteresting feature of all these ¡solated 
facts is that it has not yet been possible to establish the cause and 
effect relationship among them through a phenomonologically and mathe- 
matically defensible mechanism. In the following, the mechanism suggested 
by the numerical experiments is described in detail. 

Consider an elastically mounted circular cylinder which is at its 
lowest displacement during a given cycle. Also assume that the stationary 
cylinder value of the vortex shedding frequency is slightly greater than 
the natural frequency of the cylinder. As soon as the cylinder begins 
to move upwards, the lower attached-vortex sheet lags behind the cylinder 
because of its inertia. In other words, the wake axis passing through 
the instantaneous center of the cylinder rotates clockwise (see Fig. (49) 
for t = 189). In the meantime, the velocity relative to the cylinder 
rotates clockwise an angle dictated by the velocity of the uniform flow 
and the instantaneous velocity of the cylinder. The rotation of the 
velocity vector coupled with the rotation of the wake axis moves the 
upper separation point further downstream and the lower separation point 
further upstream relative to that which would have existed on a stationary 
cylinder. The rate of circulation at the upper separation point increases 
rapidly and decreases at the lower separation point. This gives rise to 
a larger circulation per unit time than that for the stationary cylinder 
for the vortex towards which the cylinder is moving. Furthermore, the 
vortex spiral which is attached to the upper separation point does not 
move out of the cylinder's way because of ¡ts inertia and consequently 
comes closer to the approaching cylinder. These two facts, namely the 
increase of circulation and the decrease of the relative distance between 


The cylinder and the upper vortex spiral, confirm the observation made by 
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Fig. 49. Wake configuration for self-excited transverse oscillations. 
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Fig. 49 (con't.). Wake configuration for self-excited transverse oscilla- 
ons. 
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Davies [78] that "The growing vortex on vibrating cylinders seems to roll 
up more quickly and with a strength about 35 percent larger than that for 
a stationary cylinder." The increase of the vortex strength predicted 

by the current model will be discussed later. The phenomenon is, however, 
relatively more involved than the simplified introductory remarks would 
suggest. There is a third factor which enhances the strengths of the 
attached vortex sheets. This is directly related to the difference between 
the period of oscillation of the cylinder, Ta = BÆR = Es and the time 
required for the shedding of a vortex from a stationary cylinder, ¡.e., 

e - T . Each growing vortex about an oscillating cylinder continues to 
be fed acia! not only at a higher rate but also for a longer time 
period. This further reinforces the strength of each growing vortex. 

IT is now necessary to establish the link between the enhanced 
vortex strengths and the locking-on of the vortex shedding frequency to 
The frequency of cylinder oscillation. The results of the stationary 
cylinder have shown that the larger the strength of the vortices in the 
near wake the smaller their velocity relative to the cylinder. That is, 
strong vorfices with their equally strong images tend to conserve their 
moment of-circulation. Consequently, the vortices are shed at larger 
time intervals. This, in turn, leads fo smaller Strouha! numbers. The 
decrease of 5, or the increase of le brings T closer to Tas This ten- 


dency is reinforced every succeeding cycle until Ik, =T = Je or until 


© 
the effects which cause the changes in amplitude, period, lift force, 
and the phase angle disappear. At such time, a stable periodic state is 
reached. 


The foregoing explanation dealt only with the kinematic aspects 


of the fluid motion. In the final analysis, it is the balance between 
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the work done by the fluid on the cylinder and the work done by the dissi- 
pative forces (internal friction) which determine the stability of the 
oscillations. If the net work becomes negative, then the oscillations 
decay and vice versa. The parameter which establishes the relationship 
between the friction mechanism and the given state of oscillation is the 
phase angle between the exciting force and displacement. IT is given by 
[79] 

22 0/0 

tan o = — (50) 

I- (w/w) 
As w/w approaches unity from w/w < 1, d increases rapidly, the rate of 
increase being a strong function of Z. For a cylinder with very smal | 
internal damping, Y ¡increases from say 20 degrees to d = 150 degrees for 
a very smal! change in w/w. 

The relationship between the dynamics and the kinematics of the 
motion for the three states of oscillation and hence the reason for the 
sustained oscillations in the region of synchronization may be explained 
as follows. In the case of synchronization (see Fig. 50) the maximum 
velocity of the cylinder leads the maximum lift and the lift reaches zero 
at about the Time the cylinder reaches its peak amplitude. In other 
words, the lift force opposes the cylinder motion for a short time period 
as the amplitude of the cylinder begins to decrease from its maximum 
displacement. If for.any accidental reason the amplitude of the cylinder 
in one direction were to be increased while the cylinder was undergoing 
synchronized oscillations, the increased strength of the vortex growing 
ahead of the cylinder motion increases the period, U Mas, vin Turn, 
decreases the phase angle between the force and the displacement which 


then increases the phase angle between the force and the velocity. 
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Fig. 50. Lift coefficient and cylinder velocity versus time. 
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Consequently, the energy transferred from the fluid to the cylinder de- 
creases below that dissipated by the material damping. These sequences 
of events decrease the amplitude to the level where the energy extracted 
from the fluid over a cycle is just balanced by that dissipated by 
internal triction. 

The role played by all of the terms of the governing equations of 
motion and the energy transfer mechanism may be quantified as follows. 


Re-writing Eq. (44) 


“es | TS 
O O ee 

n + Te n = 2T So ac 
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n + (44 repeat) 


and multiplying with n and integrating over a cycle, one has 


t+T t+T t+T t+T 


ame S TS 
ne O 2 OZ . ee : 
/ indt + a / n dt er ) | nndt = 2055. | Cyndt (51) 


+ í i t 





which reduces to 


++T t+T 
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In the case of a perfect balance between the energy input [right hand side 
of Eq. (52) ] and the energy dissipated (left hand side), the above equality 


holds true. For an unstable system where the equality does not hold true, 
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Thus, E = 0.0 corresponds to steady state with constant amplitude of os- 
Silat ion; E > 0.0 corresponds to amplification; and finally, E < 0.0 
corresponds to damping of oscillations. The variation of n, n, n and ai 
in Eq. (44) as a function of time is shown in Figs. 5la-c. Figure Jla 
corresponds to the case where E > 0.0; Fig. 5lb to E = 0.0; and Fig. Sic 
to E < 0.0. Clearly, the phase angle between the lift force and the 
velocity is smaller for the case of E > 0.0 than for E = 0.0. In fact, 
one can state that the said phase angle ¡increases from a plus to a minus 
value. Only in the case of Fig. lb is the phase angle such that the 
energy imbalance is for all intents and purposes zero. Furthermore, the 
Stability of the system to small perturbations in amplitude is such that 
a small increase in amplitude leads to a decrease in phase between lift 
and displacement and thus to an increase in phase between |i ft and 
velocity. This, as noted earlier, leads to a decrease in amplitude. The 
reverse is true for a slight decrease in amplitude from a steady state 
oscillation. 
5. Parametric Study of Synchronization 

Numerical experiments have been performed for the ¢ and a, combi- 
nations given in Table |! in addition to that already presented. The 
results for each set are presented in Figs. 52 through 54. All of the 
parameters shown in these figures display similar variations with w Or 
U Furthermore, they are quite similar, at least in general form, to 
those experimentally obtained by Feng (see Fig. 47). A closer look at 
the experimental and calculated values show some differences, part of 
which may be attributable to the limitations of the numerical model and 
parr to the limitations of the experimentation. In other words, all of 
the differences to be pointed out are not entirely consequences of the 


numerical model. 
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Fig. 5la. Variation of n, n, n and CL versus time for E>0.0. 
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Fig. 52. Numerical results for self-excited transverse 
oscillations with [=0.04 and a =0.025. 
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Fig. 55. Numerical results for self-excited transverse 
oscillations with G=0.08 and a =0.05. 
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Fig. 24. Numerical results for self-excited transverse 
oscillations with z=0.04 and a 70.0125. 
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A comparison of Figs. 47 and 52 shows that the numerical model 
predicts a peak lift amplitude at a U_ value smaller than ur. In Feng's 
experiments, the lift and oscillation amplitudes reached their peak at 
nearly identical U. values. The reasons for this are not entirely clear. 
It may be conjectured that it is a consequence of the extremely low damp- 
ing encountered in Feng's experiments (tf = 0.001). In fact, mechanical 
Systems with very small damping are much more finely tuned and more pre- 
cariously balanced than those with larger internal damping. IT is there- 
fore natural to expect That The more stable systems experimented with in 
The present numerical calculations should exhibit larger differences in 
U values at which they reach their extremals. 

The ratio of The cylinder oscillation frequency to the natural 
frequency, a remains practically constant ata value slightly under 
unity in E e reported by Feng and others. In The numerical 
results, w/w gradually increases from about 0.95 to 0.99 and then 
rapidly increases toward W, A The differences between the measuremenis 
and calculations may be A the difficulty of modelling some secondary 
features of flow which adjust themselves to render w/w practically 
constant. 

The amplitudes of the lift coefficient and cylinder response cannot 
be compared with the experimental data because of the large differences 
in damping and mass parameters. However, it is interesting to note that 
the ratio of the maximum lift coefficient for the oscillating cylinders 
to the lift coefficient for the stationary cylinder (standard run) varies 
from 1.3 to 1,6 (see Figs. 48 and 52-54). This lift amplification is in 


conformity with most of the data reported in the literature [75]. The 


cylinder response shows a dramatic drop for a small change in W Or Un 
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reminiscent of the hysteritic drop in Feng's data. It should, however, be 
emphasized that the drop in n in the numerical calculations has nothing to 
do with the hysteresis resulting from the change in damping of the cylin- 
der. The sudden drop is a consequence of the changes in behavior of 

flow external to the cylinder. One further distinction must be made between 
the experimentally obtained sudden drop and that obtained numerically. In 
the experiments, the velocity of the flow (or U is increased gradually 
through the U value at which the sudden change occurs. In The numerical 
calculations, U or Wo is kept constant at each point. Consequently, what 
is observed is not a sudden change in time but a change in amplitude as 

ae pune IOn Of Wo 

The cylinder response for w > wi does not remain at a constant 
amplitude. It is because of this fact that part of the response curves 
corresponding to that region are shown by broken lines in Figs. 48 and 
52-54. A sample plot of n for Uae wi is shown in Fig. 55. 

In discussing the material damping and mass parameters, another 
Parameter denoted by Sa = c/a, has been introduced. Among others, Vickery 
and Watkins [80] have attempted to show that Ny (ny at perfect synchroni- 
zation) is uniquely determined by the value of a response or stability 
parameter, S Sarpkaya [81] has shown That MM may be fairly accurately 


represented by 
Z 
x 
Ci (wr) 


Kurz O 
o De (54) 


2 
Y0.06 + Se 


where ai is the lift coefficient for the stationary cylinder. The values 
O 
of í and a, were judiciously selected to arrive at identical values of So 


for three sets of calculations and twice the SE value for one set of 
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calculations. According to the hypothesis advanced by Vickery and Watkins, 
Ty should be uniquely determinable by Se: The results obtained through 

use of the model are tabulated below. The results in Table III show that 
cylinder response at perfect synchronization is not very sensitive to the 


individual variations of the damping and mass ratios. 


TABLE 111 

Ms la 2121 
SET| ac `G 
r cocoa oar dr: 0.52 0.54 
1! | 0.04} 0.0250 | 1.6 0.58 0.58 
1111 0.08} 0.050 | 1.6 0.52 0.61 
Iv | 0.04} 0.0125 | 3.2 0.24 0.25 





An exact relationship between San and other parameters may be 


G 
established through the use of the energy balance at perfect synchroniza- 


tion. Re-writing Eq. (52) 


Eu t+T 


UNER 
O A ER 2 ° 
i J Mo = 21 So a. / C, nat (52 repeat) 


O 
t t 





and expressing £ and Ci as 


= Xx i * = Xx f * * 
en SIMI Ch Ca SIMS REDE) 


where 4% is the phase angle between the force and the displacement at 


perfect synchronization one has 


2 





TAS CE, *sinp* 
oo t+T W NM 


or 

+ Xci + 
a En 
M WS. 


Equation (56) shows that nn IS nor a unique funcTion of Se: However, 
since 9%, C, 5, and w*/w% do not vary significantly, ný varies approxi- 
mately as a A comparison of Eqns. (54) and (56) shows that Eq. (56) 


requires the use of parameters which must be determined either numerical ly 


or experimentally whereas Eq. (54) requires only a knowledge of Ch and 
O 


W Or Ur 

The results obtained through the use of Eq. (56) for the four 
sets of input parameters are also shown in Table Ill. Clearly the direct 
and indirect predictions of ni are in excellent agreement, as would be 
expected. 

The noe relationship given by Eq. (56) is expected to hold true 
for values of W. other than we provided W is confined to the region of 
synchronization. Any small deviafions between The predicfions of 

_ mw SE, sing ie 


Am W SE 


from those directly calculated with the model are simply due to small 


modulations in the amplitude of oscillations. Numerical calculations which 
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can easily be performed through the use of Figs. 48 and 52-54 have shown 
that the predictions of Eq.(56a) are within a few percent of those obtained 
directly. 

The transverse oscillations of a cylinder in uniform flow increase 
the affective projected area of the cylinder and hence the bluffness of 
the body. Consequently, hydroelastic oscillations of a cyiinder, such as 
the strumming of a cable, not only cause alternating stress, fatigue, 
and noise but also increased drag. 

The numerical model calculated two important characteristics of 
the fluctuating drag force: the mean drag and the amplitude of oscillation. 
Figure 56 shows the mean drag coefficient, Ch, and the amplitude of the 


fluctuating component, AC, 


Fig. 56 is a data point obtained experimentally [81 ].for the corresponding 


„das ec UNCION Of W. Or UL: Also shown in 


a and Ny values. 

Figure 57 is a representative plot of the instantaneous value of 
Cy as a function of time. If shows, among other things, that the drag 
force oscillates at a frequency twice that of the lift force (for lift 
i= a0 and for drag T = 3). 

The separation points on an oscillating cylinder oscillate with 
a frequency equal to that of the cylinder. Representative results pre- 
sented in Fig. 58 show that the mean position of the separation angle 
remains practically constant at a value of 9. = 77 degrees. The amplitude 
of the fluctuations follow a pattern very similar to that of Ne me 
Stagnation point undergoes similar oscillations. The results have shown 
that the instantaneous position of the stagnation point nearly coincides 


with that which would have resulted with a uniform flow at the instantaneous 


total velocity. In other words, the small oscillations of the stagnation 
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— Sarpkaya [81] 


Fig. 56. Drag coefficient versus a for 2=0,02, a 70.0125. 
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point relative to the uniform flow are obscured by the larger oscillations 
resulting from the oscillations of the cylinder or of the ambient flow 


relative to the cylinder. 
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VII. CONCLUSIONS 


The results presented herein warrant the following conclusions: 
I. A discrete-vortex model based on the rediscretization of shear layers 
and the use of appropriate vorticity cancellation mechanisms can produce 
results which are in conformity with those obtained experimental ly; 
2. The numerical model allows experimentation with the variation of any 
one of the independent parameters and permits one to understand the 
governing physical mechanisms; 
5. It has been shown that two major flow parameters, namely the drag co- 
efficient and the Strouhal number, are least sensitive to changes in the 
characteristics of the wake. The Strouhal number is intrinsic to the 
shape of the bluff body and its magnitude cannot be varied more than a 
few percent by varying the strength of the vortices by physical means. 
Experiments over the past 50 years are in conformity with these conclu- 
sions; 
4. The lift coefficient is the one most sensitive to the changes of 
flow in The near wake. This finding of the numerical model and The 
difficulty of obtaining consistent laboratory data for lift are in con- 
formity with each other. The numerical model explains the physical 
mechanism as to why lift should be a strong function of the evolution 
of the near wake; 
5. The discrete vortex model, as all other numerical methods attempting 
to model a turbulent region of a flow and its interaction with a laminar 
region, is in need of a hypothesis to properly simulate the circulation 


dissipation and the vorticity diffusion caused by turbulence. 


153 


o + 





Notwithstanding this fact, a very simple decay mechanism is capable of 
Simulating the said consequences of turbulence; 

6. The predictive powers of the numerical model have been put to test 
through its use in the study of self-excited oscillations of elastically 
mounted circular cylinders. The results were found to be in conformity 
with practically all observations and measurements made in recent years 
and led to the understanding of the mechanisms governing the synchroni- 


zation phenomenon. 
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